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Abstract Accreting compact objects are crucial to understand several important astrophysical phenomena 
such as Type Ia supernovae, gravitational waves, or X-ray and y-ray bursts. In addition, they are natural lab- 
oratories to infer fundamental properties of stars, to investigate high-energy phenomena and accretion pro- 
cesses, to test theories of stellar and binary evolution, to explore interactions between high-density plasma 
and very strong magnetic fields, to examine the interplay between binary evolution and dynamical interac- 
tions (in the case they belong to dense star clusters), and they can even be used as a probe for the assembling 
process of galaxies over cosmic time-scales. Despite the fundamental importance of accreting compact ob- 
jects for astrophysics and recent progress with the comprehension of these fascinating objects, we still do 
not fully understand how they form and evolve. In this chapter, we will review the current theoretical status 
of our knowledge on these objects, and will discuss standing problems and potential solutions to them. 
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1 Introduction 


Compact objects are remnants of stellar evolution and they are mainly divided into three categories, 
which strongly depends on the mass and metallicity of their main-sequence (MS) star progenitor. Low-mass 
stars ($ 2 Mo) and intermediate-mass stars (~ 2 — 8 Mo) end their lives as white dwarfs (WDs), while more 
massive stars evolve into either neutron stars (NSs) or black holes (BHs). Understanding stellar remnants 
represents a crucial part of understanding stellar evolution and stellar structure. Compact objects, especially 
as part of a binary, offer the potential to study physical phenomena in the most extreme environments and 
produce some of the most fascinating events in the Universe. 

Compact objects as members of binaries are frequently found to accrete material from their companion 
stars. Accreting compact objects are then binaries in which a compact object accretes matter from a com- 
panion, which can be a star, another compact object or even a planet in some cases. They are excellent 
laboratories for studying accretion disk physics and the generation/evolution of extremely strong magnetic 
fields. These objects are also among the most luminous high-energy sources in the Universe, emitting de- 
tectable radiation from X-rays to y-rays, thereby permitting detailed studies of plasma physics and the related 
emission processes under extreme conditions. 

Not only the secular evolution of accreting compact objects is intrinsically related to mass transfer, also 
their formation is in many cases associated with at least one episode of mass transfer between the stars 
composing the progenitor binary. Unlike single stars, the stars within a binary may interact and in turn may 
be affected by one another depending on their nature and the orbital separation. Mass transfer in a binary 
through the inner Lagrangian point occurs when at least one star fills its Roche lobe, which is triggered by 
either nuclear evolution of this star, for example the expansion of a giant star, or by orbital shrinkage due 
to orbital angular momentum loss (AML). In the former case, mass transfer can be maintained because the 
star expands at least at the same rate as its Roche lobe radius. The latter occurs in all binaries because two 
masses orbiting their common centre of mass always loose angular momentum due to gravitational wave 
radiation. Depending on the nature of the mass loosing star, orbital AML mechanisms that are stronger than 
gravitational wave radiation can be present. The loss of angular momentum reduces the separation between 
the two stars and therefore the radius of the companion can continuously fill its Roche lobe. 

Regardless of what is causing mass transfer, the Roche lobe filling star has to lose mass because it tries 
to keep itself inside its Roche lobe radius, which implies further mass loss than that originally expected 
from single stellar evolution. Mass transfer in binaries affects the evolution of the stellar components. For 
example, the mass growth of the helium core of a giant star can be halted before helium burning occurs if the 
giant star is loosing its entire envelope to the accreting companion and this is why most helium core WDs 
are found to be members of binaries. In other cases, the accreting star can rejuvenate by becoming more 
massive than it was at the onset of mass transfer, which alters its natural evolution. Such a mass gain could 
allow, for instance, a star that would otherwise evolve into a WD (or NS) to evolve into a NS (or BH), or 
more generally a more massive compact object than in the case of single star evolution. Another example 
corresponds to the non-expected chemical abundances in some stars such as barium stars, which can only be 
explained if their companion were mass loosing stars that managed to synthesise those elements. 

In this chapter we will discuss the state-of-the-art of our understanding of the formation and evolution 
of accreting compact objects. We start with a very brief review of the different types of accreting compact 
objects and the main processes that produce their X-ray emission. Their formation, the current state, and 
the future evolution depend fundamentally on the mechanisms that are driving mass transfer in binaries. We 
therefore subsequently review the present understanding of the formation and evolution of these systems 
with a discussion of the modes and the stability of mass transfer, before we discuss their possible formation 
channels and their secular evolution. We will then end this review with a discussion of the major current 
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problems in our understanding of the formation and evolution of accreting compact objects and briefly 
describe potential avenues towards solving these issues. 


2 The Accreting Compact Object Zoo 


Cataclysmic variables (CVs) are close binaries composed of WDs accreting from low-mass Roche lobe 
filling MS stars and are definitely the best-studied and most numerous type of interacting compact binaries. 
While several thousand CVs are known, the most important samples currently available are volume limited 
[1, 2]. These representative samples provide key observational constraints on CV formation and evolution 
theories but at the same time we have to admit that they are still rather small and to some degree subject 
to low number statistics. The majority of CVs are non-magnetic [e.g. 1] and most non-magnetic CVs are 
dwarf novae, which exhibit repetitive outbursts due to thermal-viscous instabilities in their accretion disks 
[e.g. 3, 4, 5]. All types of CVs are X-ray sources. In the case of non-magnetic CVs, the X-ray emission is 
produced in the transition region between the accretion disk and the WD, the boundary layer, from where 
approximately one half of the accretion energy is expected to be released [6, 7]. In magnetic CVs, matter 
is accreted onto the WD in a field-channeled accretion flow, starting at the region where the magnetic field 
captures the mass flow from the secondary star, and extending down to the surface of the WD. This accretion 
flow is supersonic when it reaches regions close to the surface of the white dwarf and a shock is formed. The 
matter in the post-shock region, which is the region between the shock and the WD surface, is compressed 
and heated to temperatures up to a few tens of keV. Bremsstrahlung from this post-shock region is responsible 
for most of the X-ray emission of magnetic CVs [8, 9]. When a certain amount of hydrogen-rich matter 
is accumulated on the surface of an accreting WD, hydrogen shell burning is ignited. Non-magnetic and 
magnetic CVs experience unstable hydrogen burning which produces nova eruptions [e.g. 10] and because 
of the associated mass loss, the WD masses in CVs are supposed to remain approximately constant over a 
nova cycle [e.g. 11]. 

While mass transfer in CVs is driven by orbital AML, binaries consisting of a WD and a slightly more 
massive Roche lobe filling MS star experience thermal time-scale mass transfer. In such systems, accre- 
tion rates exceeding values of ~ 1077 — 1076 Mo yr7!, depending on the WD mass, can be reached and the 
hydrogen-rich material accreted by the WD can stably burn on the surface of the WD [e.g. 12]. Nova erup- 
tions and the corresponding erosion of the WD mass are avoided in theses systems and therefore they are con- 
sidered to be potential Type Ia supernova progenitors [13]. Hydrogen burning on the surface of a WD is pro- 
ducing effective temperatures somewhat below 100 eV, producing X-ray emission in the 0.01—2 keV range. 
Stably hydrogen burning WDs in close binaries are therefore called super-soft X-ray binaries (SSXBs). The 
first two SSXBs were discovered in the eighties [14] and about a decade later one additional system was 
found [15]. Today the discovery of more than 100 systems has been reported in external galaxies, the Mag- 
ellanic clouds, and the Milky Way, but this number and the total emission in X-rays seems to be too low to 
account for the full rate of Type Ia supernovae [16, 17]. The relatively low number of SSXBs and the corre- 
sponding low soft X-ray emission is potentially related to the fact that circumstellar matter can suppress the 
soft X-ray flux [18]. 

A binary where a compact object is accreting material at a relatively high rate (greater than several 
107? Mo yr7!) from a low-mass evolved red giant is called a symbiotic star (SySt), in case the compact 
object is a WD [19], or a symbiotic X-ray binary (SyXB), in case the accretor is a NS [20]. In most cases 
(~ 80%), the donor is a normal evolved giant, typically M-type red giant, but in some cases it can be a 
G-/K-type yellow giant. Their near-infrared emission shows the presence of a stellar photosphere, and for 
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this reason these systems are called S-type SySt [21, 22, 23, 24]. In the remaining systems, the donor is a 
Mira-type variable and their near-infrared emission is consistent with being originated from a combination 
of a reddened Mira-type variable and a warm optically thick dust shell, which makes these systems to be 
called D-type SySts [25]. In SySts, the accretion rates are typically large enough to cause stable nuclear 
burning on the surface of the WD, which leads to the emission of soft X-rays similar to those of SSXBs. In 
most SySts and SyXBs, the mass transfer from the red giant to the compact object is due to gravitationally 
focused wind accretion and only in some cases due to atmospheric Roche lobe overflow. SySts are important 
for understanding mass transfer in wide binaries and the formation of astrophysical jets. Most importantly, 
as SSXBs, SySts are promising candidates for being Type Ia supernova progenitors [26, 27, 28]. 

AM CVn stars are compact binaries with orbital periods in the range of ~ 5— 65 minutes [29]. These 
binaries are so close that a MS star cannot be part of the system as it would just not fit. In AM CVns, a WD 
accretes helium-dominated material from a Roche lobe filling companion, usually a lower mass WD or a 
semi-degenerate helium-rich star. Mass transfer in AM CVns is driven by orbital AML through gravitational 
wave radiation which is expected to be detectable by the Laser Interferometer Space Antenna (LISA) mission 
[30, 31]. Observationally, AM CVns are characterized by their hydrogen deficient optical spectra and their 
blue colour. More than 50 AM CVns are known by the time of writing [e.g. 29] and their number is constantly 
increasing thanks to dedicated spectroscopic and photometric transit surveys [32, 33]. X-ray emission can be 
triggered in two different ways in AM CVns and the details depend to some degree on the formation channel 
and the nature of the donor. In case the donor is a WD, the onset of mass transfer occurs as a direct impact 
during which a small region of the accreting WD is heated and emits a fraction of the released accretion 
luminosity at high energy wavelength [34]. In the presence of an accretion disk around the accreting WD, 
X-ray emission is produced in the boundary layer in full analogy to the X-ray emission of non-magnetic 
CVs. 

Low-mass X-ray binaries (LMXBs) consist of a NS or a BH that is accreting from a Roche lobe filling 
low-mass star [35]. As in CVs, mass transfer in LMXBs is driven by orbital AML due to magnetic brak- 
ing and/or gravitational wave radiation. In analogy to dwarf novae, LMXBs also show outbursts that are 
most likely triggered by thermal-viscous instabilities in their accretion disks [36]. An important difference 
between LMXBs and CVs is that in the former irradiation of the accretion disk is more important and may 
fully suppress disk instabilities. In CVs, disk irradiation is dominating most likely only in post-novae sys- 
tems [37, 38]. LMXBs also show thermonuclear X-ray bursts originating from unstable hydrogen and helium 
burning on the surface of the accreting NS [39], similarly to the nova eruptions occurring in accreting WDs. 
A typical LMXB emits almost all of its radiation in X-rays and these systems are therefore among the bright- 
est X-ray sources. Since the seventies, an over-abundance of LMXBs in globular cluster in comparison to the 
Milky Way field has been observationally established, which naturally suggests that dynamical interactions, 
which are very common in the dense central parts of globular clusters, facilitate their formation [e.g. 40]. 

Ultra-compact X-ray binaries (UCXBs) are binaries with periods shorter than about an hour. UCXBs are 
very similar to AM CVns with the main difference being that the accreting compact object is either a NS or 
a BH. As in the case of AM CVns, such short orbital periods exclude MS stars to form part of the binary and 
therefore the donors are either WDs or semi-degenerate helium-rich stars [41]. UCXBs are mainly observed 
at X-rays and some of these objects show thermonuclear X-ray bursts [42]. As in the case of LMXBs, there 
is currently a clear consensus that the relative number of UCXBs is much larger in globular cluster than in 
the Milky Way field [40], and dynamics therefore seem to provide an important formation channel also for 
these systems. 

If a compact object is accreting from a high-mass star (= 10 Mo), the system is known as a high-mass 
X-ray binary (HMXB) [43]. As in the cases above, the potential energy of the accreted material in HMXBs is 
converted into X-ray emission. Although HMXBs are therefore strong X-ray sources, they are also bright at 
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optical wavelengths given that their optical emission is dominated by the massive companion. In some cases, 
the donor is an O- or early B-type supergiant that is filling or close to filling its Roche lobe and these systems 
are consequently called supergiant HMXBs (sg-HMXBs) [e.g. 44]. Despite hosting a supergiant, the orbital 
periods of some sg-HMXBs are shorter than ~ 10 days, and they are persistent X-ray sources, powered by 
the capture of matter from the strong stellar wind of the supergiant or, in a few cases, through Roche lobe 
overflow [45]. In contrast, supergiant fast X-ray transients (SFXTs) have orbital periods typically longer that 
~ 10 days and exhibit peculiar characteristics compared to the remaining sg-HMXBs, such as bright short 
time-scale flares. A second sub-type of HMXBs are B-emission (Be-HMXBs), which were discovered in 
the 70ies [46]. In virtually all cases these systems contain accreting NSs, but there is one confirmed Be- 
HMXB hosting a BH [47], and a few Be-HMXBs have been proposed to contain WDs [e.g. 48, 49]. The 
donors in these systems are rapidly rotating B stars on or close to the MS, and unlike the donors in sg- 
HMXBs, given their typical orbital period (~ 10—400 days), the donors in Be-HMXBs are deep inside their 
Roche lobes [e.g. 50]. Be-HMXBs can be in slightly eccentric orbits [51], or highly eccentric orbits, and 
their X-ray emission is originated when the compact object accretes from the Be star circumstellar decretion 
disk either during periastron passages, resulting in orbital-phase-dependent X-ray outbursts lasting several 
days, or due to enhanced activity of the Be star, which leads to X-ray outbursts lasting several weeks not 
correlated with any particular orbital phase [52]. As a final sub-class, we mention Wolf-Rayet HMXBs 
(WR-HMXBs), in which the donor is a massive helium star. Only few such systems are known and they 
are likely descendants from sg-HMXBs and Be-HMXBs [53]. Among HMXBs, Be-HMXBs correspond to 
the most numerous class (~ 60% [54]), followed by ss-HMXBs (~ 32% [54]). This is a direct consequence 
of selection effects favouring the discovery of Be-HMXBs in comparison with sg-HMXBs. As previously 
mentioned, Be-HMXBs are discovered during outbursts, which faint persistent ss-HMXBs would not have. 
For this reason, it is very likely that the contribution of sg-HMXBs to the intrinsic population of HMXBs is 
much higher. 


3 Modes of mass transfer 


The Roche lobe is the last equipotential surface in a binary that keeps all the mass attached to one star. 
Any expansion of a star beyond its Roche lobe or any orbital shrinkage of the Roche lobe due to orbital 
AML may cause a star to fill its Roche lobe thereby forcing mass to flow over to the other star. Determining 
the rate at which mass is transferred through the inner Lagrangian point in binaries represents in general a 
difficult hydro-dynamical problem. One of the most frequently used and reasonable prescriptions properly 
takes into account the finite scale height of the stellar atmosphere of the Roche lobe filling donor. 

For stars with extended atmosphere Roche lobe overflow can occur through the inner Lagrangian point 
even when the donor is slightly under-filling its Roche lobe, i.e. Rg < Rap, where Rg is the donor radius, RRL 
is the Roche lobe radius of the donor. Taking this into account, the mass transfer rate Ma can be estimated 
as [55] 

(1) 


: Ret — R 
Ma = -Mo exo(- a 7 ; 


Ap / yq) 
where Hp is the pressure scale height at the photosphere of the donor. 


The quantity y(q) is a function of the mass ratio q = Ma/Ma, where Ma is the donor mass and M, is the 
accretor mass, given by 
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7 0.954 + 0.025 logjoq — 0.038 (logjoq), if 0.04<q<1 iy 
MD = 10.954 + 0.039 logigg + 0.114 (logigq”’, if 1<q<20 
Finally, Mo is the mass transfer rate when the donor exactly fills its Roche lobe given by 
. Qn Rè \/keTea y” 
Mo = Fop | 7 | z 3 (3) 
V2.71828 / | GMa J \ mp Hph 


where G is the gravitational constant, kg is the Boltzmann constant, mp is the proton mass, Mg is the donor 
mass, Teg is the donor effective temperature, uph and Ppp are the mean molecular weight and density at its 
photosphere, respectively and F(q) = 1.23 + 0.5log,,q, valid for 0.5 < q < 10. Outside the ranges of validity, 
F(q) and y(q) are usually evaluated using the value of q at the edge of their respective ranges [56]. 

In case the donor is over-filling its Roche lobe, i.e. Ra > Rp, the mass transfer rate can be expressed as 
[57] 


. : RÈ PRL 1/2 2 r kgT 1/2 
Ma = -Mo - 27 F(q) ri —]| qP, (4) 
GMa Poh I,+1 Mp Ll 


where T is the temperature, J”; is the first adiabatic exponent, and Pph and Pry are respectively the pressures 
at the photosphere and at the radius corresponding to Rpt. 

We shall mention that, in case donor is strongly under-filling its Roche lobe, the mass transfer rates pro- 
vided by Equation | are negligible. However, mass can still be transferred from the donor, if the stellar winds 
from the donor are sufficiently strong and dense and a part of the this wind is accreted by the companion. 
In what follows, we discuss in more detail these modes of mass transfer, starting with the stability of mass 
transfer when the donor is filling its Roche lobe, i.e. Ra ~ RRL. 


4 Stability of mass transfer through Roche lobe filling 


To understand the evolution of binaries containing a Roche lobe filling star it is of fundamental importance 
to understand the reaction of the donor and the Roche geometry to mass loss. Despite the fact that the 
constraints on a Roche lobe filling star are modified due to mass loss and the star consequently can no longer 
be considered as evolving in isolation, under certain conditions the star may adjust its internal structure 
in response to mass loss because most of its mass is concentrated towards its centre. Whether hydrostatic 
and/or thermal equilibrium can be maintained in a Roche lobe filling star and whether mass transfer is stable 
or unstable depends entirely on the response of the donor radius and its Roche lobe radius to mass loss on 
the corresponding time-scales [58]. 


4.1 Dynamical time-scale mass transfer 
Changes of the radius of the donor star as a response to variations of its mass and through time can be 


expressed as 


dln Ra = 


ðln Ra ( ðln Ra 


ar * (amm) Me 8 
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where ¢ is the time. While the first term on the right side is connected with stellar evolution and relaxation 
processes connected with thermal readjustment, the second term corresponds to the adiabatic response of 
the donor to mass loss, i.e. the hydrodynamic response of the donor radius to variations in its mass. 

The changes in the donor Roche lobe radius as a response to variations in mass, time and angular mo- 
mentum can be written as 


ðln RRL 


InR 
dln RRL = + k ZARI 


ðln Ma 


PT ) dln Ma. (6) 
The first term on the right side represents the changes in the Roche lobe radius due to orbital AML in the 
absence of mass transfer, and the second term is associated with the response of the Roche lobe radius to 
mass loss. 

Whether mass transfer from a Roche lobe filling star is stable or not can be evaluated by comparing two 
logarithmic derivatives associated with the donor. These derivatives are that of the donor radius with respect 
to its mass, for a fixed initial profile in composition and specific entropy, and that of its Roche lobe radius 
with respect to its mass. The former is usually called adiabatic radius-mass exponent ćaq, defined as 


ðln Ra 
ad E 7—7] o (7) 
oln Ma ad 
while the latter is termed the Roche lobe radius-mass exponent Cp, defined as 
dln Rep 
= —__, 8 
éRL = ayn Me (8) 


Combining Equations 5 and 6, and using the definitions 7 and 8, we can express the change in the donor 
radius relative to its Roche lobe radius as 


ain( x = ( oan z ae dt + (Ga-Sei)din Ma . (9) 

The condition for the donor to remain within its Roche lobe in hydrostatic equilibrium during an episode 
of nearly stationary mass transfer is obviously dln(Ra/R) ~ 0. Given that mass transfer is occurring, the 
difference in the time derivatives in Equation 9 must be positive, because the donor radius needs to change 
more strongly than its Roche lobe radius to guarantee that the donor keeps filling its Roche lobe. In addition, 
as the donor is losing mass, dln Mg is necessarily negative. Therefore, the difference between the adiabatic 
and Roche lobe radius-mass exponents needs to be positive to keep the donor in hydrostatic equilibrium. 
In other words, by comparing these two derivatives we are able to determine whether the response to mass 
transfer of the donor inner layers can occur adiabatically, i.e. whether the entropy profile in the stellar interior 
remains unchanged despite the mass loss. If, on the other hand, Zaq < Gay, the mass transfer is unstable and 
occurs on the dynamical time-scale, because the mass transfer rate is only limited by the sonic expansion of 
the envelope through the inner Lagrangian point. 

Another way of reaching the same conclusion is by analysing perturbations of the mass transfer rate [e.g. 
55]. During nearly steady mass transfer, i.e. dIn(Ra/RL) ~ 0, the well-known mass transfer rate equation as 
a function of ćaq and ¿RL can be derived from Equation 9 as 


Ma ea OlnRey 
ad — SRL 


(10) 


i --( at at 
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Considering a system in an unperturbed phase of stationary mass transfer, i.e. Ma ~ 0, we can evaluate 
the dynamical stability of mass transfer after a sudden non-null perturbation 5Mq. After the perturbation, 
Ma= Mao +6Mq, where Mao is the unperturbed mass transfer rate, and Ma œ 6MaMa (Cad — SRL). Given that 
Mg < 0, the mass transfer will be dynamically stable only if the response to a sudden increase in |Ma] is a 
decrease in |Ma]. This happens only if ća — č&RL > 0. Otherwise, the sudden increase in |Ma] will result in 
further increase of |Ma], which leads then to dynamically unstable mass transfer. 

Even though realistic and detailed binary calculations must be performed to investigate the stability of 
mass transfer [e.g. 59, 60, 61, 62], polytropes can shed some light on this process in a very simple way 
and provide a reasonable first-order approach to the problem in some cases [e.g. 63, 64]. More specifically, 
composite polytropes, which are supposed to reproduce rather well the structure of fully convective low- 
mass MS stars, provide ćaq in very good agreement with those computed from detailed stellar models [62]. 
On the other hand, condensed polytropes, which are approximations of red giants, are unable to reproduce 
the non-ideal gas effect and inefficient convection occurring in detailed stellar models, which implies that 
condensed polytropes likely underestimated ¢aq of red giants [62]. If simple stellar models do not provide 
realistic estimates, the adiabatic radius-mass exponent needs to be calculated by solving the differential 
equations for stellar structure. 

Once Zaq has been obtained, it needs to be compared to ¢RL in order to evaluate whether mass transfer is 
dynamically stable or not. The Roche lobe radius-mass exponent Zp, can be expressed as 


Olna ðln(RgL/a) Olng 
= — 11 
SRL OlnMg i ( Olng ln My!’ i) 
where a the semimajor axis. The three terms in Equation 11 can be written as [65, 64, 60] 
l 2g = i= 
ðlna _ TE |! +qU-fhm)| 4 (12) 
OlnMg q+1 q+1 
ðln(RrL/a) 2 (q'/3\[1.2g"/3 +1/0 +q!) (13) 
ðlng 3 3 }| 0.6g2/3 +n + g!/3) J’ 
and 
Olng 
=1 1- — 14 
din Mg +( Onl Bm) 4. (14) 


where 0) < œm < 1 corresponds to the fraction of mass lost from the vicinity of the donor in the form of fast 
isotropic winds, and 0 < Bm < 1 corresponds to the fraction of mass lost from the vicinity of the accretor 
due to isotropic re-emission. Given that the fraction of mass leaving the donor that is effectively accreted 
is 1 — @m — Bm, the fully conservative case corresponds to @m = Bm = 0, while ami + mi = 1 represents 
the fully non-conservative case. In a more general situation, ¿r depends on several other processes, which 
are related to orbital AML and mass loss from the system [64], such as emission of gravitational waves, 
magnetic wind braking, tidal interactions, and torques exerted by circumbinary disks. The most important 
dependencies of Zp with respect to the stability of mass transfer are those on the mass ratio as well as on 
the fraction of mass lost from the binary. 

The relation between ¢aq and ¿RL for a binary composed of a red giant donor with initial mass 1.2 Mo and 
a companion of mass 1.1 Mo is illustrated in the left-hand panel of Figure 1 (taken from a binary model [60]). 
As mentioned before, the limit for dynamical stability strongly depends on mı. In particular, the greater Bm, 
the smaller Zp_. This means that the greater mı, the better the chances for the red giant donor to remain 
in hydrostatic equilibrium. If wind mass loss is included in the simulations, r, decreases even further for 
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more evolved red giants (i.e. red giants with larges core masses), for a given Bm. For instance, when mass 
loss through winds is not considered, fully conservative mass transfer is dynamically stable only for core 
masses > 0.63 Mo, but this limit becomes = 0.46 Mo, if wind mass loss is included. 


T T T T T 7 T £ T T T T T T 
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Fig. 1 Adiabatic (¢,4) and Roche lobe (gL) mass—radius exponents for a red giant donor of initial mass 1.2 Mo with a com- 
panion of mass 1.1 Mo, as a function of the red giant core mass (Ma,c), taken from a binary model (left-hand panel [60]), and 
the critical mass ratio separating dynamically stable from unstable mass transfer for a donor of mass 5 Mo, as function of its 
radius (Rq) during its evolution, taken from three different binary models (right-hand panel [59, 61, 62]). In the computation of 
RL in the left-hand panel, for six different fractions of mass lost due to isotropic re-emission from the vicinity of the accretor 
(Bm1), mass loss through stellar winds are considered (red lines) or not (blue horizontal lines). In the right-hand panel, wind 
mass loss is not considered in either binary model and fully conservative mass transfer is assumed. It is clear from the left-hand 
panel that non-conservative (mı > 0) mass transfer is more easily dynamically stable in comparison with the fully conservative 
case (Bm = 0), regardless of whether wind mass loss is taken into account or not. This happens because érg becomes smaller 
as Bm increases, which allows mass transfer to be dynamically stable for smaller core mass, i.e. the more mass loss through 
isotropic re-emission is allowed, the more extended the range of stability. The inclusion of wind mass loss enhances this effect, 
by decreasing even further ¢p_, for core masses greater than ~ 0.3 Mo, extending in turn even further the range of stability. 
The right-hand panel, on the other hand, clearly illustrates how the assumptions made in binary models can drastically change 
the boundary separating dynamically stable from unstable mass transfer. For instance, the solid critical mass ratio strongly 
increases during the AGB evolution, which is not seen in the other two models. In addition, the dashed critical mass ratio is on 
average much greater than the other two models. 


Even though the fundamental principles related to dynamical stability of mass transfer are well under- 
stood, outcomes of binary models are strongly dependent on the physical assumptions. As mentioned earlier, 
for a given binary at the onset of mass transfer, ér, strongly depends on the mass ratio and on the assumed 
mass loss and orbital AML prescriptions. In addition, ¢,q depends entirely on the donor mass and its evolu- 
tionary status. Thus, for a given binary, it is possible to derive a critical mass ratio separating dynamically 
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stable from unstable mass transfer, as a function of the evolutionary status of the donor. This critical mass 
ratio is plotted as a function of the donor radius on the right-hand side of Figure 1, for a donor with mass 
5 Mo, taken from three different binary models [59, 61, 62]. Even though all binary models show the initial 
decrease in the critical mass ratio as the donor evolves along the first giant branch (FGB), it is clear from 
the figure that these models provide qualitatively and quantitative different results. In particular, one of them 
delivers on average a much higher critical mass ratio, in comparison with the other two. In addition, while in 
two models the critical mass ratio slowly decreases as the donor climbs the asymptotic giant branch (AGB), 
the other model predicts a rapid increase of the critical mass ratio. The main differences in the three works 
are connected with the mass transfer rate equation (whether it is limited by the local sound speed or not), 
the radius used to quantify ¢aq (surface or inner) and the Lagrangian radius adopted to quantify ¢py (inner or 
outer). 

Therefore, despite significant progress during the last decade, the limiting mass ratio for the onset of 
dynamically unstable mass transfer remains quite uncertain. This has significant implications for our general 
understanding of the formation of close binaries hosting compact objects, as dynamically unstable mass 
transfer is considered a crucial ingredient in most formation scenarios proposed for these objects. If mass 
transfer is dynamically unstable, there are mainly two outcomes. If the donor is either an MS star or a 
degenerate object, then most likely the binary will merge on the dynamical time-scale of the donor, with a 
non-null fraction of the total mass being eventually lost in such an energetic event. On the other hand, if 
the donor has evolved beyond the MS, common-envelope (CE) evolution will be triggered, resulting in a 
strong decrease of the orbital period. In fact, CE evolution triggered by dynamically unstable mass transfer 
represents the standard scenario for the formation of close binaries containing a compact object. Without 
a proper understanding of the conditions for dynamically unstable mass transfer, it will be impossible to 
reliably predict the relative numbers of objects as important as SN Ia progenitors and binaries that are main 
gravitational wave sources. 


4.2 Thermal time-scale mass transfer 


If the Roche lobe filling donor in a binary can maintain hydrostatic equilibrium, the natural follow- 
up question is whether it must sacrifice its own thermal equilibrium or not to restrain its surface within 
its Roche lobe. The thermal equilibrium of a star, i.e. the balance between the nuclear energy production 
and the atmospheric radiation losses, is disturbed as a response to mass loss. However, as in the case of 
hydrostatic equilibrium, the donor may be able to re-establish thermal equilibrium, if mass transfer proceeds 
on a time-scale slower than its thermal time-scale, i.e. if the donor interior has enough time to relax to 
thermal equilibrium. Otherwise, mass transfer will be thermally unstable, i.e. it will proceed on the thermal 
time-scale of the donor. In order to verify whether this happens or not, we use the thermal equilibrium 
radius-mass exponent eq, defined by 


Z dln Ra 
“4 din Maln ` 


In order to calculate ¢eq, the mass-radius relation for stars in thermal equilibrium must be known. If mass 
transfer is dynamically stable, it will also be thermally stable if ¢eq > rx. In this case, mass transfer needs 
to be driven by an external process and proceeds on a time-scale different from the dynamical and thermal 
time-scales. On the other hand, if faq > ¿RL > Seq, Mass transfer will be driven by thermal readjustment and 
proceed on the donor thermal time-scale. 


(15) 
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4.3 Nuclear or orbital angular momentum loss time-scale mass transfer 


If mass transfer is dynamically and thermally stable, continuous mass transfer can still be generated by 
either stellar or binary evolutionary processes. In this case, mass transfer is not self-stimulated as in the 
previous cases, since mass transfer does not proceed by the virtue of the inability of the donor to remain in 
either dynamical or thermal equilibrium. Dynamically and thermally stable mass transfer is driven by either 
nuclear expansion of the Roche lobe filling star or orbital AML. In the former case, the donor keeps filling 
its expanding Roche lobe, while in the latter case, the Roche lobe and the donor shrink as a response to mass 
loss. In case mass transfer is driven by nuclear expansion of the donor, mass transfer proceeds at a rate such 
that Ra ~ RgL. On the other hand, if mass transfer is driven by orbital AML, then the mass transfer rate is 
entirely dependent on the strength of orbital AML. The main orbital AML mechanisms driving the evolution 
of binaries are gravitational wave radiation (GR) [e.g. 66] and magnetic wind braking (MB) [e.g. 67, 68]. 

Using the weak-field approximation of general relativity, it is possible to derive the orbital angular mo- 
mentum changes due to emission of gravitational waves from two point masses as [69] 
1+ ge | = 


Jor Mı M2(M, + M2) 


4 


(1-2) 


where a is the semimajor axis in Ro, Mı and M3 are the binary component masses in Mo, e is the eccentricity, 
and Jom is the orbital angular momentum given by 


ia = a mys |í z ) ee (17) 


(16) 


orb a 


= -8.315x107!° | 


Mı +M2 Porb 


where Por» is the orbital period in years. Emission of gravitational waves acts all the time in every binary. 
However, from the strong dependence on the orbital separation, only for very close binaries it will have any 
impact on the orbital angular momentum evolution. Of course, how close a given binary must be so that GR 
is non-negligible depends on the component masses and on the eccentricity. 

While Equation 16 is well-established and gravitational wave radiation is well understood, the same 
cannot unfortunately not be said about MB. The main physical processes that drive MB are quite easy to 
understand but as usual the devil is in the details. It is clear that a highly ionized wind leaving the star is 
forced to corotate with the magnetic field of the star out to the Alfvén radius, thereby taking away more 
angular momentum than the specific angular momentum of the star which consequently spins down. If the 
star is tidally locked to the orbit, i.e. if the star is synchronized with the orbital motion, the mass loss through 
winds removes significant amounts of orbital angular momentum, even if the wind mass loss rate is rela- 
tively small. The torques produced by the magnetized winds depend on the strength of the magnetic field, 
the Alfvén radius, the wind mass loss rate, and on the stiffness of the magnetic field at large distances from 
the star. However, the details of these dependencies are very uncertain and neither observationally nor the- 
oretically well constrained. Consequently, in the last decades several attempts to formulate empirical laws 
or recipes that are applicable to particular cases have been developed and a large number of MB prescrip- 
tions is available on the market [70]. Unfortunately, as we will see, for a given fixed star mass, these MB 
prescriptions provide dramatically different orbital AML rates spanning easily several orders of magnitude. 

In what follows, we briefly mention some important MB prescriptions, but discuss later in more detail 
their impact on the evolution of accreting compact objects. A MB prescription that is widely used for stars 
in open star clusters and based on several simplifying assumptions for the winds, magnetic field geometry, 
and flow acceleration profile, has been developed by S.D. Kawaler [71] 
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where ag depends on the scaling between the secondary magnetic field and spin, n depends on the secondary 
magnetic field geometry, Mwina < 0 is the star mass loss rate through winds, Q, R and M are the star spin, 
radius and mass, respectively. Even though the parameter Kw depends on the ratio of the wind speed to the 
escape speed at the Alfvén radius and how the star magnetic field scales with the star radius and spin, it 
is usually approximated by ~ 2.035 x 10°3(24.93”)(6.63 x 10°)#”/3. With ax = 1 and n = 1.5, Equation 18 
describes the standard Skumanich law [72], as orbital AML due to MB becomes proportional to 2°. 

On the other hand, in the context of accreting compact objects driven by orbital AML, usually the so- 
called RVJ prescription [73] is adopted, which is given by 


l M\/Ro\¥[ Q2\(R\® 
Juprvi = -38x100 . \(2) (2) (=) dyn cm, (19) 


where 0 < ymp < 4 is a free parameter that largely affects how Jp varies during the binary evolution. With 
ymB = 4, Equation 19 describes the standard Skumanich law. 

Despite the fact that the RVJ prescription can explain several observations of accreting objects, it has 
serious problems in different contexts, mainly connected with the mass transfer rates it provides and the 
corresponding evolutionary time-scales during CV and LMXB evolution [e.g. 74, 75, 76, 77]. In order to 
overcome part of these difficulties, recently a new convection and rotation boosted (CARB) prescription that 
enhances orbital AML due to MB was proposed [78, 79]. The CARB prescription can be expressed as [79] 


11/3 8/3 pags 
‘ 2 7 -1/3 4/3 8/3 Q T, 2 2 Q- R 
JMB,CARB = — 3 (—Mvina) RP Qo BY Qe = Vese + K > (20) 


where vesc is the surface escape velocity, Tconv is the convective turnover time, B is the surface magnetic field 
strength, and K2 = 0.07. The main difference between the RVJ and the CARB prescription is that the former 
is the result of an empirical fitting scheme, while the latter was obtained in a self-consistent physical way, 
considering wind mass loss, the dependence of the magnetic field strength on the outer convective zone and 
the dependence of the Alfvén radius on the donor spin. 

We illustrate in Figure 2 the orbital AML rates provided by the three MB prescriptions mentioned above. 
The evolutionary sequences have been computed by us using the Modules for Experiments in Stellar Astro- 
physics (mesa) code [80, 81, 56, 82, 83], in which we evolved a Roche lobe filling star with initial mass of 
1 Mo, orbiting a point-mass companion of mass 1.4 Mo. In the simulations, we take into account stellar wind 
mass loss, adopting the Reimers prescription [84] with the scaling factor equal to 1, and assume fmi = 0.8 
and @m = 0. For MS masses < 0.4 Mo, the CARB prescription clearly provides higher orbital AML rates 
than the RVJ and the Kawaler prescriptions. For MS masses greater than that, the CARB, RVJ and Kawaler 
prescriptions provide comparable orbital AML rates if nx = 1.5 and ywp = 0. Finally, GR is typically much 
weaker than any of these prescriptions, expect for the Kawaler recipe with n < 1. 

In addition to GR and MB, there are other sources of orbital AML acting on interacting binaries which 
arise as a consequence of mass transfer. For this reason this type of orbital AML is called consequen- 
tial AML. Among the several mechanisms proposed for consequential AML, the most important are a fast 
isotropic wind from the vicinity of the donor, isotropic re-emission from the vicinity of the accretor [85], 
a circumbinary disk [86] and winds from the accretion disk [87]. Consequential AML is by definition an 


Contents 15 


1038 ee |, rc ff °° - yg |  « « «fo. 


CARB (subgiant) 
si CARB (MS) ao saa E use ins satents keene 
1037 a a, 4 
Kawaler ae 
(n = 1.5) 
10° 4 
= 
d 
= ! 
S 10% f 4 
= l 
K: i STRESS zer E a Tr a a a a e tele lt 
1034} i Leeri a J 
i ae 
a S E S E O E rs 
1033 = | Mm fF 4 J er | 
ce ne eee e a a oo, 
ae cece a Ce a Gravitational Wave Radiation 
1075 i 4 
D i S G a a es ee ee 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
M (Mo ) 


Fig. 2 Comparison of the orbital AML rates (—J) due to GR and due to some MB prescriptions, during the evolution of a low- 
mass Roche lobe filling star, as a function of its mass (M). All tracks have been computed by us with the mesa code, assuming 
solar metallicity, wind mass loss (Reimers parameter equals to 1), Bm = 0.8, @m = 0, a point-mass accretor with initial mass 
of 1.4 Mo, and the initial mass of the donor is 1.0 Mo. Each track represents the evolution from the onset of Roche lobe filling 
to the moment when the star becomes fully convective (when MB is supposed to cease or at least significantly decrease). The 
dot-dashed red lines corresponds to the CARB prescription [79], for which we chose two different initial orbital periods, namely 
12 hours to illustrate the case of an unevolved MS star (helium core mass equals to zero), and 4 days to illustrate the subgiant 
case (helium core mass ~ 0.1 Mo). The solid blue lines shows the RVJ prescription [73], for three different values of ymp, and 
in all these cases we chose an initial orbital period of 12 hours, corresponding in turn to unevolved MS stars. The Kawaler [71] 
prescription is marked by the long-dashed green lines, for three different values of n, assuming in all three sequences ax = 1, for 
which we also chose an initial orbital period of 12 hours, corresponding to the unevolved MS star case. Finally, the short-dashed 
black lines correspond to the strength of GR in each of the evolutionary sequences, which can be easily associated with them 
by comparing the mass at which the MB becomes inefficient, i.e. the lowest mass in each track. It is clear that the CARB recipe 
provides on average higher orbital AML rates than the RVJ and the Kawaler prescriptions, the latter for n < 1.5, especially for 
masses < 0.4 Mo. On the other hand, the CARB, RVJ and Kawaler prescriptions provide comparable orbital AML rates for 
larger masses, if nx = 1.5 and ymp = 0. In addition, the Kawaler prescription is much more sensitive to the chosen parameters, 
in comparison with the RVJ prescription, as it provides much more different orbital AML rates, for small variations of nx. 
Moreover, except for the Kawaler prescription with n < 1, GR in all evolutionary sequences contribute only negligibly to the 
total orbital AML. Finally, regarding the type of the star in the CARB tracks, we can see that the orbital AML rate is initially 
higher for a more evolved star, but decreases as the star mass decreases, resulting in a weaker MB at masses < 0.4 Mo, in 
comparison with an unevolved MS star. 


orbital AML mechanism that is mass-transfer dependent, i.e. it is null in the absence of mass transfer. It is 
usually written in a simple parameterized form as 
JCAML Ma 


=v—, (21) 
Jorb Ma 
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where v depends on the binary properties, including perhaps the mass transfer rate itself, and is the parameter 
characterizing the strength of consequential AML. Its form depends entirely on the mechanism driving 
orbital AML. 

In the case of a circumbinary disk, v is given by [86] 


ven = Ym (1 +q), (22) 


where 6m is the fraction of the mass lost from the binary components that feeds the circumbinary disk, and 
the parameter y takes into account the characteristics of the circumbinary disk as a function of time, and 
depends on the inner edge of the circumbinary disk, the binary semimajor axis, the viscous time-scale at the 
inner edge of the circumbinary disk and the mass of the circumbinary disk. 

In the case of mass loss from the binary, v may correspond to the specific angular momentum of the 
accretor (isotropic re-emission) or the donor (fast isotropic wind), given by [85] 


q 1 
Vml,a = Smi 1+4 > Vml,d = ml T+q : (23) 


Fast isotropic winds and isotropic re-emission may also generate additional orbital AML due to friction 
arising from the interaction of the ejected material and the accretor/donor. This frictional AML then strongly 
depends on the velocity of the ejected material and on the orbital velocity of the accretor/donor as well as on 
the geometrical cross-section of the accretor/donor [88, 89, 90, 91]. 

We have already illustrated in Figure 2 that MB typically provides orbital AML rates several orders 
of magnitude higher than those expected from GR. The rates predicted by consequential AML may be 
higher or lower in comparison with MB, depending on the form of the two types of orbital AML. For 
instance, frictional AML due to isotropic winds and/or isotropic re-emission is typically negligible, unless 
the expansion velocity of the ejected material is very small [92]. The same also occurs for consequential 
AML due to mass loss (isotropic winds or isotropic re-emission), which typically provides rates smaller 
than those related to either MB or GR alone. On the other hand, consequential AML due to a circumbinary 
disk can be the dominant source of orbital AML depending on the parameters of the disk and the rate at 
which it is fed [93, 94]. 

So far we have briefly described the current understanding of the different processes that can drive mass 
transfer in binaries. This included the stability of a given binary against mass transfer on the dynamical and 
thermal time-scales of the donor as well as mass transfer driven by the nuclear evolution of the donor and 
orbital AML. The largest uncertainties in current models are probably related to the conditions separating 
dynamically stable from dynamically unstable mass transfer, the strength and parameter dependencies of 
orbital AML through magnetic braking, and the processes that drive consequential AML. On the basis of 
the above-described knowledge of mass transfer in binaries, formation channels for the different types of 
accreting compact binaries have been developed in the last decades which we will review in the following 
section. 


5 Formation Channels 


Most accreting compact objects form from the evolution of zero-age MS—MS binaries. Some may form 
from stellar mergers in hierarchical systems, i.e. from triples, quadruples and, eventually, higher order sys- 
tems. Dynamical interactions, which are common in high-density environments such as the central parts of 
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dense star clusters, could also lead to the formation of a sufficiently important fraction of accreting com- 
pact objects, especially for those originating from high-mass stars. While for some specific systems, and 
some particular types of accreting compact objects, the triple and/or the dynamical scenario is required to 
explain their properties, in most cases the contribution from these two scenarios is expected to be much 
smaller than that of binary evolution. This is simply because of the lower frequency of triples and higher or- 
der hierarchical systems in comparison to binaries [e.g. 95] and because strong dynamical interactions only 
occur frequently in dense star clusters. In what follows, we therefore mainly focus on accreting compact 
objects formed though binary evolution, but add a few comments on additional formation channels through 
dynamical interactions at the end of this section. 

The current small orbital separations of close binaries in which a compact object accretes from a non- 
degenerate low-mass companion, i.e. CVs and LMXBs, imply that these systems are currently much smaller 
than the progenitor of the compact object was as a giant. Consequently, the orbit of the binary must have 
been significantly reduced during the formation of the compact object, most likely through CE evolution 
generated by dynamically unstable mass transfer [96, 97]. 

The current orbital periods of most wide binaries in which a WD accretes from the winds of an evolved 
low-mass red giant, i.e. most S-type SySts, are significantly longer than those of post-CE binaries but short 
enough that mass transfer must have occurred when the WD was formed. Therefore, in these cases, the most 
likely scenario is that the orbit expanded during the formation of the WD, which happens when the mass 
transfer is dynamically stable and most likely non-conservative [98]. 

SyXBs are similar to SySts, in the sense that they are also wide-orbit binaries in which the compact object 
accretes from the wind of an evolved low-mass red giant. However, their compact objects are NSs, which 
implies that stable non-conservative mass transfer cannot be the formation channel, because the required 
mass ratio of the zero-age MS-MS binary implies that the non-degenerate companions should be high-mass 
stars, which is not the case. A possible formation scenario for these systems could be the formation of a long- 
period post-CE binary, hosting a helium star paired with a low-mass star, followed by either electron-capture 
or core-collapse supernova leading to the NS [99, 20]. 

Regarding compact objects accreting from high-mass non-degenerate stars, i.e. HMXBs, those in wide 
orbits (e.g. Be-HMXBs) are most likely formed through an episode of dynamically stable non-conservative 
mass transfer, while those having orbital periods of only a few days (e.g. sg-HMXBs) seem to have formed 
through CE evolution [100]. On the other hand, NSs/BHs accreting from the strong and dense winds of 
Wolf-Rayet stars, i.e. WR-HMXBs, require two episodes of mass transfer, one to form the accreting compact 
object and the other one to form the Wolf-Rayet star. Even though the first episode of mass transfer is most 
likely dynamically stable non-conservative mass transfer [101], there is evidence that in some WR-HMXBs 
the second episode of mass transfer was CE evolution, while others can be more easily explained when the 
second episode was dynamically stable non-conservative mass transfer [53]. 

Finally, in those very close binaries in which the compact object is accreting from a Roche lobe filling WD 
or semi-degenerate helium-rich star, i.e. AM CVns and UCXBs, two episodes of mass transfer are needed, 
and at least one of these episodes must have been CE evolution. More specifically, these systems may form 
either through two simply CE phases, or through dynamically stable non-conservative mass transfer followed 
by CE evolution, or finally through CE evolution followed by dynamically stable non-conservative mass 
transfer [102, 103]. 

In what follows, we will in more detail address the main formation channels of accreting compact objects, 
which are illustrated in Figure 3. 
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Fig. 3 Main formation channels leading to accreting compact objects. Starting from a zero-age MS—MS binary, when the more 
massive MS star evolves and fill its Roche lobe, mass transfer stars. Mass transfer is dynamically stable if the adiabatic radius- 
mass exponent (Z,q) is greater than the Roche lobe radius-mass exponent (Za), and unstable otherwise, leading to CE evolution. 
During dynamically stable mass transfer, the orbital period increases or decreases as a response to mass transfer, and the orbital 
period at the end of mass transfer depends on the binary properties at the onset of Roche lobe filling and on the fraction of 
the transferred mass that is lost (@m) + mı + mı) during the event. On the other hand, if CE evolution happens, drag forces 
remove orbital energy, which is used with a certain efficiency (@cg) to eject the envelope, potentially together with a fraction 
of other sources of energy, such as recombination energy (Œ@rec) and/or thermal energy (ain). The outcome of CE evolution is a 
binary with a much shorter orbital period. Compact objects in close orbits accreting from a Roche lobe filling non-degenerate 
low-mass companion (i.e. CVs and LMXBs) are believed to form during one episode of CE evolution. NSs accreting from 
the winds of low-mass red giants in wide orbits (i.e. SyXBs) are also expected to form in an episode of CE evolution. On the 
other hand, WDs in wide orbits accreting from the winds of low-mass red giants (i.e. SySt) are expected to form through an 
episode of dynamically stable non-conservative mass transfer. Regarding NSs or BHs accreting from high-mass companions 
(i.e. HMXBs), they can be formed through dynamically stable non-conservative mass transfer or CE evolution. Finally, compact 
objects accreting from Roche lobe filling WDs or semi-degenerate helium-rich stars (i.e. AM CVns and UCXBs) are formed 
through two episodes of mass transfer, at least one of them being CE evolution. 
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5.1 Common-envelope evolution 


When the more massive star of a zero-age MS—MS binary evolves off the MS and fills its Roche lobe, 
it will give rise to dynamically unstable mass transfer, if Cag < aL. In this case, the mass loss time-scale 
becomes so short that the donor cannot remain within its Roche lobe, which leads to the formation of a 
CE around the dense giant core and the MS star. Even though a CE event in many situations is triggered 
by dynamically unstable mass transfer, there are also other mechanisms able to do that, such as the Darwin 
instability [e.g. 104], which occurs if the future donor angular momentum at synchronism exceeds one-third 
of that of the orbit. The Darwin instability corresponds to a tidal instability. It occurs when not enough 
orbital angular momentum is extracted from the orbit to maintain the star synchronized as it evolves. Tidal 
forces will spin up the star by removing orbital angular momentum, resulting in a binary with smaller orbital 
separation and in turn even less orbital angular momentum. This implies that the star will need even more 
orbital angular momentum to stay synchronized with the orbital motion. This leads to a runaway process of 
orbital decay and ultimately a CE event. 

Because of drag forces during CE evolution, the MS stars and the core spiral in towards their common 
centre of mass and the CE is expelled. In this process, a large fraction of the initial orbital energy and orbital 
angular momentum is lost with the envelope which leads to a post-CE binary with an orbital period orders 
of magnitude shorter than that of the initial MS—MS binary. In case dynamically unstable mass transfer 
was generated from a massive star to its companion, the post-CE binary hosts a helium star, which is the 
progenitor of the compact object, and the binary separation may increase due to a natal kick during the 
collapse of the helium star. 

Despite a lot of effort and (computational) resources that are spent towards a better understanding of CE 
evolution, we still struggle to identify and describe some of the main physical processes involved in this 
crucial evolutionary phase. The main reason for this failure lies in the fact that modelling CE evolution is 
very demanding from a computational point of view, since it involves a complex mix of physical processes 
operating over a huge range of scales [105]. In the absence of reliable simulations that cover the full com- 
plexity of CE evolution, a very simplistic approach is frequently used to overcome this problem and to make 
predictions for populations of post-CE binary characteristics. These simple equations based on energy and 
angular momentum conservation provide useful constraints on the outcome of CE evolution. 


5.1.1 The energy budget of common-envelope evolution 


In the so-called energy formalism, the outcome of CE evolution is usually approximated by the balance 
between the change in the orbital energy and the envelope binding energy. This energy conservation equation 
is parameterized with a parameter 0 < acg < 1, which corresponds to the efficiency with which orbital energy 
is used to eject the envelope. In other words, œcg represents the fraction of the difference in orbital energy 
(before and after the CE phase) that unbinds the envelope, i.e. 


GMa M, a (24) 


Evina = ace AEom = - ec | —>— - 
2af 2a; 
where Find is the envelope binding energy, Erp is the orbital energy, Ma, is the core mass of the donor, a; 
is semimajor axis at onset of CE evolution, and ap is the semimajor axis after CE ejection. 


The binding energy is usually approximated by 
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where 4 is the binding energy parameter, which depends on the structure of the donor [106, 107]. 
In case thermal and/or recombination energy stored inside the envelope are assumed to help in ejecting 
the envelope, the the binding energy equation is usually written as 


Ma Gm Ma Ma 
Ebina = -f — dm + an f En(m) dm + ane f Erec (m) dm (26) 
Mac r(m) Mac Mac 


Ebina = — (25) 


where &p and £rec are the specific thermal and recombination energies and @ and Qyec are the fractions of 
the thermal and recombination energies that are assumed to contribute to unbinding the envelope. Thermal 
energy is the thermodynamic internal energy of the CE, which is constrained by the Virial theorem, and can 
be released and partially used to help in the CE ejection. Recombination energy is stored in ionized helium 
or hydrogen within the CE. These ions can recombine in the expanding layers of the gaseous CE, once it has 
sufficiently cooled and expanded, thereby releasing energy [108] that can contribute to unbinding the CE. 
The binding energy parameter 2 can be computed by equating Equations 25 and 26. When œm > 0 and/or 
Qrec > 0, additional sources of energy contribute to the ejection of the CE, which implies that less orbital 
energy is needed to unbind the envelope. 

Deriving constraints on the values of @cg, Qh, and @rec has been a focus of close compact binary research 
in the last decades. Unfortunately, several slightly different CE prescriptions are available which needs to be 
taken into account when comparing the results of different studies. In what follows we briefly describe the 
constraints obtained so far for binaries hosting WDs, for which the progenitors are intermediate-mass MS 
stars, and for those hosting NSs or BHs, which originate from the evolution of high-mass MS stars. 

Regarding low-/intermediate-mass stars, the best way for constraining the CE evolution is by using sam- 
ples of detached WD-MS post-CE binaries because of the large number of known systems and because it is 
relatively easy to characterize them [109]. Reconstruction of the evolution of observationally characterized 
post-CE binaries consisting of WDs with M-type companions suggests that CE ejection is rather inefficient, 
with œcg ~ 0.2 —0.3 [110]. Comparisons between predictions from binary population synthesis with obser- 
vations reached similar conclusions for this type of post-CE binary [1 11, 112, 113]. Recently, observational 
surveys have become capable of identifying close WD binaries with MS companions more massive than 
M-type MS stars (spectral types A, F, G and K) and early results indicate the existence of a population of 
post-CE binaries that can be reproduced with the same small CE efficiency [1 14, 115]. However, for some 
of these systems, a small fraction of recombination energy seems to be required to contribute to CE ejection 
[116]. 

Close binaries hosting hot subdwarf B stars (i.e. helium-core-burning stars with very thin hydrogen-rich 
envelopes), mostly paired with WDs but some also with M-type MS stars, are also reasonably good targets 
for estimating the efficiency of CE ejection, as CE evolution is certainly involved in their formation. Contrary 
to WD-MS post-CE binaries, results from binary population synthesis of hot subdwarf B stars indicate that 
the observations are better reproduced with acg ~ 0.75 and a large fraction (~ 0.75) of the internal energy 
(thermal plus recombination) contributing to expelling the envelope [1 17, | 18] which indicates that both CE 
ejection and the conversion of thermal and recombination energies into kinetic energy are efficient processes. 
However, given that the sample of WD-MS binaries is cleaner in the sense that only one phase of mass 
transfer is involved, the observational constraints indicating a small value for the CE efficiency should be 
considered as being stronger. In addition, we have to be aware that the energy prescription for CE evolution 
is too simplistic and several unknown physical processes are not properly taken into account which might 
affect the constraints derived from observations of different systems on the CE efficiency. 
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With respect to progenitors of accreting compact objects hosting NSs and BHs, the situation is more 
complicated because we know less about CE evolution in this mass range [105] and the situation might be 
quite different from the previous case. The main reason for these expected differences is that CE evolution 
leading to WD-MS post-CE binaries is generated by dynamically unstable mass transfer from red giant 
donors with a convective envelope [119], while for high-mass stars with solar metallicity the onset of Roche 
lobe filling usually happens when the donor is in the Hertzsprung gap and has a radiative envelope. In 
addition, metallicity has a strong impact on the evolution of high-mass stars prior to the onset of Roche lobe 
overflow, which is independent of the effect of metallicity-dependent winds, thereby altering the evolutionary 
stage at which Roche lobe filling occurs [120]. In most cases, this causes 4 to be far lower for high-mass 
stars than for intermediate-mass stars [121], implying that in several situations, there would be apparently 
not enough energy available to eject the CE. 

For a long time, this has been believed to be the case for BH-LMXBs, which would require an abnormally 
high efficiency for the envelope ejection in binary models [e.g. 122], since it is very difficult for the binary 
hosting a Roche lobe filling high-mass star paired with a much less massive companion to eject the CE at 
the expense of only orbital energy. However, more recent calculations have shown that this is not an actual 
issue in the formation of BH-LMXBs, since they could be produced without invoking unrealistically high 
CE efficiencies [123]. The real problem arises when trying to reproduce the donor mass distribution in this 
population. Despite the fact that it is possible to reproduce individual BH-LMXBs in binary models, we 
cannot yet reproduce the properties of the entire population with just one model [123]. We will come back 
to this issue later. 

Regarding SyXBs, in the absence of other sources of energies, in population synthesis aiming at explain- 
ing their observational properties, «cg needs to be as high as 4 [20], which is unrealistically high. This result 
could suggest that there are other sources of energy are required to help to expel the envelope or that these 
X-ray binaries do not form through standard CE evolution . For instance, they could instead be preferentially 
formed from accretion-induced collapses of a massive oxygen-neon-magnesium WD [e.g. 124]. Another 
possibility, perhaps more realistic, is that the uncertainties and simplifications involved in stellar and binary 
evolution calculations, for example mixing in stellar interiors, core-collapse/electron-capture supernova and 
CE evolution, are the main reason for the problems in reproducing these systems, as well as BH-LMXBs 
[e.g. 125]. 


5.1.2 Common-envelope evolution from hydro-dynamical simulations 


Despite the above described attempts to constrain @cg, th and Q@rec, it seems unlikely that there is a 
universal recipe within the simplistic energy formalism able to account for all types of accreting compact 
objects. In the case of low-/intermediate-mass stars, evidence is growing evidence that the efficiency of 
expelling common envelopes is rather low and overall the energy formalism delivers reasonable results. In 
contrast, it is very difficult to explain the properties of accreting compact objects coming from high-mass 
stars based on the energy balance of CE evolution, regardless of the assumed values for the CE evolution 
parameters. This difficulty must to some degree arise from our inability to understand the role played by 
other sources of energies acting together with orbital energy, as well as correctly computing the CE binding 
energy in an unbiased way. 

CE evolution can be divided into several subsequent phases, which helps to better understand why this 
event is so complicated to be modeled. Most of the problems arise from the fact that each of these phases 
occurs at a particular time-scale, involving a particular physical process [e.g. 126]. The first phase is the 
loss of co-rotation, which can last hundreds of years, and corresponds to very small changes in the orbital 
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separation in comparison with the orbital period. During this phase, the accretor orbits either outside the 
future CE, or inside the CE outer layers [e.g. 127]. The second phase is called plunge-in, which corresponds 
to the rapid spiral-in with huge changes in the orbital separation. During this phase, the accretor plunges 
inside the CE and the orbital energy is transferred to the CE, driving its expansion, in case there is enough 
energy available for that. At the end of this phase, most of the CE mass is outside the orbit of accretor and the 
core of the donor [e.g. 127]. During the third phase, when the CE is well extended, a self-regulating spiral-in 
occurs on the thermal time-scale of the CE, during which energy released due to the spiral-in is transported 
to the CE surface and subsequently radiated away. This phase corresponds to small changes in the orbital 
energy [e.g. 127]. CE evolution terminates at the end of the self-regulating spiral-in phase with either the 
ejection of the CE or the accretor or the core of the donor filling its Roche lobe. The latter most likely result 
is a merger, but could also provide another pathway for CE ejection [e.g. 128]. 

We illustrate CE evolution in the left-hand panel of Figure 4, where we show the evolution of the orbital 
separation from just before the plunge-in phase until just after the self-regulating spiral-in, for two models 
taken from a study based on a three dimensional smoothed particle hydrodynamics code [127]. In these 
simulations, the authors chose a low-mass red giant donor with a mass of 1.8 Mo, a core mass of 0.318 Mo, 
and radius of 16.3 RO, and several companion masses. The figure clearly illustrates that there is a huge 
orbital shrinkage during the plunge-in phase until the onset of the self-regulating spiral-in phase, which is 
precisely what the energy formalism tries to reproduce. However, in those cases that avoid merger, the CE 
is not always effectively ejected which represents a long-standing problem in numerical simulations of CE 
evolution [e.g. 130, 131]. 

Given the above-mentioned inefficient ejection of the CE predicted in simulations, there should exist at 
least one energy source, perhaps more, in addition to orbital energy, that would be able to expel the remaining 
bound matter [e.g. 131]. It has been argued for a long time that thermal energy and mainly hydrogen and 
helium recombination energy could solve this problem [e.g. 132, 133, 134]. The impact of including thermal 
and recombination energies is illustrated in the right-hand panel of Figure 4, which shows the evolution of the 
unbound mass fraction through time, taken from a numerical simulation with a moving-mesh hydrodynamics 
code [129]. In this simulation, the donor is on the AGB undergoing helium-shell burning and no thermal 
pulses, with a mass of 0.97 Mo and a core mass of 0.545 Mo. Almost the entire CE mass is unbound when 
the stored internal (thermal plus recombination) energy is used, most of which is only released at rather late 
times, when the expanding CE sufficiently cools down. This is not the case when only the potential energy 
contributes to the binding energy, since ~ 10% of the CE mass is still inside the orbit of the donor core and 
the companion, after the plunge-in phase. No further dynamical spiral-in is expected, because of the low 
density of the remaining envelope and the co-rotation of the material, avoiding in turn any drag force. 

The example above, together with other simulations, provide support to the claim that internal (mainly 
recombination) energy might be a key ingredient in CE evolution. On the other hand, most of the available 
recombination energy could be simply radiated away [e.g. 135, 136], since the optical depth is expected 
to be low when recombination occurs, allowing radiation to more easily escape. In addition, despite the 
fact that the amount of available recombination energy depends on the core definition, it is expected to be 
smaller in high-mass stars stars than in intermediate-mass stars. Therefore, some simulations suggest that 
recombination energy is unlikely to help in the removal of the CE in this mass range [137]. 

In any event, it is widely accepted by now that only dynamical spiral-in is most likely not enough to 
completely remove the CE. In case recombination energy is not significantly helping, other sources of energy 
would be needed [e.g. 131]. These additional contributions could come from nuclear energy [128], accretion 
energy [138], dust-driven winds [139], interaction of the binary (core of the donor + companion) with a 
circumbinary disk [140], large-amplitude pulsations [141], or jets launched by the companion as it accretes 
mass from the circumbinary envelope when it is about to exit the CE from inside [131]. Even though there are 
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Fig. 4 Change in the orbital separation (left-hand panel [127]) and envelope mass (right-hand panel [1 29]) through time during 
CE evolution, taken from two different CE evolution models. In the left-hand panel, the onset of the plunge-in phase (when 
(|4|Porb)/a becomes greater than 0.1), the end of the plunge-in phase (when (|@| Pop) /a becomes smaller than 0.1), and 
the onset of the self-regulating spiral-in (when |(Eorb Port ) / Eorp | becomes smaller than 0.01), are indicated with crosses. The 
transition between the end of the plunge-in phase and the onset of the self-regulating spiral-in is not sudden in these simulations, 
although it lasts for only a few days. In the model with companion mass of 0.1 Mo, there is not enough orbital energy available 
to expel the CE and the binary is expected to merge during the self-regulating spiral-in phase. However, in the model with 
companion mass of 0.2 Mo, even though orbital energy is transferred to the CE, it is not enough to entirely eject the CE 
during the self-regulating spiral-in phase. Despite the fact that a strong reduction in the orbital separation occurs, like in the 
0.2 Mo model, in most numerical simulations, the CE is not entirely ejected, which indicates that other sources should contribute 
to the CE ejection together with orbital energy. In the right-hand panel, there is a comparison between the case without internal 
energy and the case with internal energy (thermal plus recombination), which has been widely proposed as the key source of 
energy usually missing in the simulations. It is clear in this simulation that the CE can be ejected with the help from internal 
energy. However, there are strong physical arguments suggesting that this energy would be simply radiated way, not actually 
contributing to the CE ejection. If confirmed beyond doubts, other sources of energy or transfer mechanisms must be playing a 
role and should be further explored as alternatives to recombination energy. 


several mechanism available, they have been poorly explored in hydro-dynamical simulations, and it is still 
not clear which one is dominant, if any at all, and under which conditions. This means that CE evolution is a 
subject far from being well understood and, hopefully, in the next years, further hydro-dynamical simulations 
will shed more light on this extremely important phase for the formation of accreting compact objects. 
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5.2 Dynamically stable non-conservative mass transfer 


Even though the energy formalism is successful to some extent in avoiding the messy physics involved in 
the several phases of CE evolution, it fails to reproduce some observed systems, in particular double helium 
WD binaries. Explaining these systems with two consecutive CE events requires that energy is generated 
during CE (which implies that the energy conservation law is violated), or at best, that an unknown source 
of energy, much greater than the orbital energy, is contributing to expelling the envelope (which can be 
considered unrealistic). 

To solve this problem, an alternative but similarly simplistic formalism was invented based on conser- 
vation of angular momentum, which is usually called angular momentum (or y) formalism [142, 143]. Al- 
though this new approach could somehow work mathematically, it does not provide a physical explanation 
for the required CE phases without spiral-in [e.g. 132, 105]. 

We believe that most systems for which no reasonable solutions in the framework of CE energy balance 
exist, are actually not formed through CE evolution, but rather dynamically stable non-conservative mass 
transfer. Indeed, there is strong theoretical evidence that an increase in the orbital separation can be reached 
during an episode of dynamically stable non-conservative mass transfer on a thermal and/or nuclear time- 
scale, which can potentially reproduce otherwise inexplicable systems. In what follows, we discuss in more 
detail systems that may form through this dynamically stable non-conservative mass transfer channel. 


5.2.1 Low-/intermediate-mass stars 


There are several types of binaries hosting (proto) WDs with orbital periods longer than ~ 100 days. 
In some of them, it is clear that some sort of interaction in their history must have occurred, because of 
observed chemical anomalies. For instance, barium stars, which are paired with WDs, are F-/G-type MS 
stars or G-/K-type red giants exhibiting strong absorption lines of ionized barium in their spectrum as well 
as of other s-process elements [e.g. 144, 145], which could not have been synthesized by them. In a similar 
way, carbon-enhanced metal-poor stars show abundances of s-process elements [e.g. 146, 147] that cannot 
be explained in the context of single star evolution. S-process elements are produced during the thermally- 
pulsing AGB phase [e.g. 148] and it is therefore widely accepted that these elements were transferred to 
the barium stars and carbon-enhanced metal-poor stars via stellar winds or dynamically stable Roche lobe 
overflow, when the WD progenitor was at this evolutionary stage. 

Other types of binaries having comparable orbital periods include, post-AGB binaries, with orbital periods 
in the range of ~ 100 — 2500 days [149]. In addition, some blue stragglers, which are typically found in old 
star clusters and correspond to MS stars bluer, brighter and more massive than those around the cluster 
turnoff point, are in binaries with orbital periods between ~ 100 and ~ 3000 days. While some hot subdwarf 
B stars are members of close binaries, a certain fraction of these stars are members of binaries with rather 
wide orbits. Hot subdwarf B stars are post-FGB binaries that must have formed when their progenitors 
filled their Roche lobe just before reaching the tip of the FGB. The typical orbital periods for wide binaries 
containing hot subdwarf B stars range from ~ 700 to ~ 1300 days. Moreover, most S-type SySts have orbital 
periods between ~ 200 and ~ 1500 days. In some cases, there is evidence for negligible mass accretion 
during these episodes of mass transfer, such as for hot subdwarf B stars [150], while in other cases, such as 
for barium stars and carbon-enhanced metal-poor stars, only a non-negligible amount of mass accretion can 
explain their chemical properties [e.g. 151, 152, 153]. 

All these binaries, in addition to having orbital periods usually longer than ~ 100 days, have in many cases 
non-null eccentricities. We show the orbital period versus eccentricity diagram for them in Figure 5, where 
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Fig. 5 Diagram depicting the orbital period (Porb) versus eccentricity for several types of low-/intermediate-mass post-mass- 
transfer binaries with orbital period longer than ~ 100 day, namely barium star binaries [154], blue straggler binaries [155, 156], 
post-AGB binaries [149], s-process-polluted carbon-enhanced metal-poor star (CEMP-s) binaries [146, 147], wide-orbit hot 
subdwarf B binaries [157, 158], S-type SySts [21, 159, 160, 161, 162, 163], and R Aqr, which is the only D-type SySt with radial 
velocity data suitable for orbital determination [25]. The hatched regions correspond to the location of post-CE binaries (Por S 
100 days) and of systems formed without Roche lobe filling (Porb 2 2000 days), taken from a binary model [164]. It is clear 
from the diagram that there is a significant fraction of systems in eccentric binaries with orbital periods between ~ 100 and 
~ 2000 days, which is precisely between the region occupied by post-CE binaries and the region occupied by systems formed 
without Roche lobe filling. Therefore, most of these systems apparently cannot be explained by CE evolution, nor by the lack 
of Roche lobe filling. Indeed, the overwhelming majority of these systems are most likely formed during dynamically stable 
non-conservative mass transfer, accompanied by an eccentricity pumping mechanism that would efficiently work against tidal 
circularization for some system, but not for all of them, as it must also account for those systems in circular orbits. 


we also added a hatched region indicating the location of the post-FGB and post-AGB binaries formed 
without Roche lobe overflow (taken from a binary model [164]). We can see from the figure that those 
systems with orbital periods longer than ~ 2000 days fit reasonably well within the expected region for WD 
binaries where the WD progenitor never filled its Roche lobe. The same certainly holds for most D-type 
SySts, which have long orbital periods (= 40 years), although only for one system the orbital period has 
been accurately measured. The mass exchange that is required to explain the formation of barium stars and 
s-process-polluted carbon-enhanced metal-poor stars with periods exceeding ~ 2000 days has likely been 
driven by efficient wind accretion, instead of Roche lobe overflow. 

In contrast, systems with orbital periods shorter than ~ 2000 days must have experienced Roche lobe 
overflow. Given their current orbital period, the progenitors of the WD (or hot subdwarf B stars) in these 
systems must have filled their Roche lobes at some point along their evolution. It is difficult to explain the 
periods of these systems as a result of CE evolution, since orbital periods of post-CE binaries are usually 
shorter than ~ 100 days [165, 166] and often just a few hours [109]. In addition, even though some of these 
binaries are on circular orbits, many others have eccentric orbits, which is not expected for post-CE binaries, 
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since tidal forces during CE evolution would clearly have managed to circularize their orbits. Consequently, 
Figure 5 indicates that a different mode of mass transfer has to be considered to understand these systems. 
For orbital periods longer than the typical outcome of CE evolution, but shorter than those that might avoid 
Roche lobe overflow, i.e. = 100 and < 2000 days, dynamically stable and non-conservative mass transfer 
represents a reasonable alternative. 

Indeed, several theoretical studies show that dynamically stable non-conservative mass transfer can ex- 
plain these systems [e.g. 98, 59, 167, 168, 169]. Depending on the combination of initial masses, metallicity, 
and orbital period of the low-/intermediate-mass zero-age MS—MS binaries, the more massive star could 
fill its Roche lobe while on the FGB (having or not a degenerate core) or on the AGB. The properties of 
post-stable-mass-transfer binaries then depend on the above-mentioned combination of initial masses, initial 
orbital period, and metallicity as well as on the fraction of mass and orbital angular momentum lost by the 
binary during the process. 

We show evolutionary tracks for dynamically stable non-conservative mass transfer in Figure 6. It is clear 
that orbital periods between a few hundreds and a few thousands days can be achieved during dynamically 
stable non-conservative mass transfer. Extrapolating the results shown in the figure towards smaller initial 
periods one can easily imagine that stable mass transfer may also lead to orbital periods < 100 days. 

According to Figures 5 and 6, formation through dynamically stable but non-conservative mass transfer 
therefore represents a reasonable scenario for most blue stragglers, post-AGB binaries, long-period binaries 
hosting subdwarf B stars, and a significant fraction of barium stars and s-process-polluted carbon-enhanced 
metal-poor stars. Also most S-type SySts, which correspond to ~ 80% of known SySts [21, 22, 23, 24], 
should be formed through an episode of dynamically stable non-conservative mass transfer. Their WD and 
red giant masses are mostly in the range ~ 0.4—0.6 Mo and ~ 1.2—2.2 Mo [21], which does not seem 
difficult to obtain with dynamically stable non-conservative mass transfer, especially when taking into ac- 
count a potential WD mass growth during SySt evolution. Following the episode of mass transfer, the WD 
companion will eventually evolve beyond the MS and become an evolved red giant, either close to the tip of 
the FGB, or on the AGB. In most cases, the mass transfer from the red giant to the WD is enhanced due to 
gravitationally focused wind accretion [170, 171, 172, 173, 174, 175]. If the accretion rate onto the WD is 
sufficiently high, the binary will become a SySt. 

Despite the general agreement between predictions and observed systems, there are three S-type SySts 
(T CrB [176], RS Oph [159, 177], and V3890 Sgr [178]) that have very different WD and red giant masses, 
and their red giants are nearly filling their Roche lobe. In these three systems, the compact objects are massive 
carbon-oxygen WDs with masses in the range ~ 1.2— 1.4 Mo, and their Roche lobe filling red giants have 
relatively low masses (~ 0.7 — 1.2 Mo). These characteristics exclude that these systems formed through 
stable mass transfer as the mass of the donor star (i.e. WD progenitor) in the progenitor system was far more 
massive than the accretor (i.e. red giant progenitor). It is therefore very likely that these three SySts have 
formed through CE evolution [27], and their WD/red giant masses have substantially increased/decreased 
during their evolution. Why common envelope evolution resulted in exceptionally long orbital periods in 
these cases (~ 200 — 600 days) remains unclear. However, if additional energy (e.g. recombination energy) 
is efficiently contributing to expelling the envelope such long post common envelope periods seem possible 
[179, 166]. Following the CE evolution, the WD companion becomes an evolved red giant, its extended 
atmosphere eventually fills its Roche lobe and mass transfer proceeds similarly to what is illustrated in 
Figure 6. We will discuss in more details the secular evolution of the different types of SySts later. 

Even though dynamically stable non-conservative mass transfer is a very good candidate to explain the 
orbital periods of a large fraction of systems with orbital periods ranging from ~ 100 to ~ 2000 days, the 
observed eccentricities remain to be explained. In order to have orbital periods consistent with observations 
at the end of mass transfer, the initial orbital periods cannot be too short, nor too long. This would imply 
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Fig. 6 Evolution towards a post-FGB binary with the donor mass (Ma) of the orbital period (Porb), in the left-hand panel, 
and the mass transfer rate (Ma), in the right-hand panel, during dynamically stable non-conservative mass transfer using the 
MESA code. We chose two different initial orbital periods (100 and 300 days) and seven different values for the transferred mass 
loss fraction from the vicinity of the accretor due to isotropic re-emission (mı), from 0.1 to 0.7 in steps of 0.1, which are 
indicated in the figure (the upper curve in each set of tracks in the left-hand panel corresponds to By = 0.1). The simulations 
start with a zero-age MS-MS circular binary, with stars of masses 1.50 and 1.43 Mo, with negligible rotation. In addition, 
we considered solar metallicity, changes in the rotation of the stars, tidal synchronization as well as wind mass loss (Reimers 
parameter equals to 0.5). Furthermore, except for stellar winds, we assumed that any mass loss in the form of fast isotropic wind 
from the donor is negligible (œmı = 0), and ignored mass loss and orbital AML due to a circumbinary disk (6, = 0). We can see 
some nice correlations in the left-hand panel of this figure, in which the colours indicate the evolution of the accretor mass (Ma) 
along each track. In particular, the smaller Bm1, the longer Porb at the end of mass transfer, which occurs because the increase 
in the orbital period due to mass transfer is much stronger than the increase due to mass loss from the binary. In addition, as 
expected, the smaller mı, the larger the fraction of mass effectively accreted by accretor, and in turn, the greater M, at the 
end of the mass transfer episode. Finally, the longer the initial Porp, the longer Por» at the end of mass transfer. The right-hand 
panel of the figure illustrates the typical behaviour in such episodes of dynamically stable mass transfer [e.g. 60], namely the 
transition from a thermally unstable to a thermally stable phase of mass transfer, the later being driven then by the nuclear 
nuclear evolution. The mass transfer starts thermally unstable because the equilibrium radius-mass exponent (¢eq) is smaller 
than the Roche lobe radius-mass exponent (Zr), which means that it is driven by the thermal readjustment of the donor, and 
proceeds initially on the donor thermal time-scale. However, as the donor loses mass and changes its structure, eq eventually 
becomes larger than Zp, and mass transfer proceeds on the donor nuclear time-scale. We can also see a nice correlation in this 
panel, the smaller Bm), the greater the peak of Ma during thermally unstable mass transfer. This happens because the donor is 
increasingly driven out of thermal equilibrium, as mass transfer tends to be more conservative. On the other hand, the smaller 
Pm, the lower the mass transfer rate during thermally stable mass transfer, i.e. the slower the mass transfer, which has an impact 
on the post-FGB mass. In particular, the smaller Bm, the larger Ma at the end of mass transfer, which can also be seen in the 
left-hand panel and happens because the overall evolution is slower when mass transfer tends to be more conservative, allowing 
in turn the core mass to grow more as Bm decreases. 
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then that, even for an initially eccentric orbit, tidal interaction should be strong enough to circularize the orbit 
before the onset of mass transfer. However, it is clear from Figure 5 that many systems have eccentric orbits 
at the end of mass transfer. This leads then to three possibilities, either we still do not fully understand the 
tidal interaction phenomenon, specially the circularization time-scales, or there is an eccentricity pumping 
mechanism operating during mass transfer, or both effects are acting together. 

Concerning the first possibility, i.e. tidal forces, recent observational efforts on ellipsoidal red giant bi- 
naries in the Large Magellanic Cloud suggest that tidal interaction should be much weaker than expected 
from theoretical considerations [180]. Among the 81 binaries that were investigated, 20% have eccentric 
orbits, which should not exist according to our understanding on tidal dissipation in convective stars. In or- 
der to explain such systems, the circularization efficiency needs to be two orders of magnitude smaller than 
predicted. 

Regarding the possibility that the large eccentricities are generated by some kind of eccentricity pump- 
ing, we note that there are currently several mechanisms proposed that could at least help to explain the 
non-circular orbits of the post-stable-mass-transfer binaries. First, the eccentricity can increase if the mass 
loss rate from the donor is higher near periastron than near apastron. The required enhanced mass loss oc- 
curs during periastron passages can be caused by either wind accretion or Roche lobe overflow or even a 
combination of both, in which there is a smooth transition between wind accretion and Roche lobe filling 
[181, 182, 157]. 

Second, in the circumbinary disk scenario, eccentricity pumping occurs due to the resonances in the 
interaction of the binary with a circumbinary disk [183, 157]. Since circumbinary disks are most likely 
created only from mass lost from the outer Lagrangian point, the eccentricity pumping due to circumbinary 
disks should actually work together with the phase-dependent Roche lobe overflow scenario, by enhancing 
even further the eccentricity. A third possibility is that the proto WD receives a natal kick after the mass 
transfer episode, similarly to natal kicks in the formation of NSs, although not that strong. 

The final possibility is the triple scenario, in which dynamical and stellar evolution with a tertiary com- 
panion would be responsible for creating eccentric binaries hosting compact objects, via chaotic orbital 
evolution of the stars, which can trigger close encounters, collisions, exchanges between the stellar compo- 
nents [e.g. 184]. Even though all these scenarios are promising, and provide reasonable results for particular 
systems, they all fail to some extent to explain the overall properties of the systems and further investigations 
are needed [e.g. 157, 185, 186, 187]. 

There is an important question to be answered regarding the picture we described above. If the accretion 
efficiency during dynamically stable non-conservative mass transfer is high, what happens with the accretor 
spin? This issue has been addressed in binary models [188], which suggest that only ~ 0.05 Mo of accreted 
material is usually enough to drive the accretor to critical rotation. This is a major problem to explain barium 
stars and s-process-polluted carbon-enhanced metal-poor stars unless there exists a mechanism allowing the 
accretor to transfer its spin angular momentum back to the orbit with the help of tidal torques from the donor 
[188]. However, it is not yet clear what kind of mechanism could do this job. 


5.2.2 High-mass stars 


After discussing dynamically stable non-conservative mass transfer among low-/intermediate-mass stars, 
we now turn to the implication of this mass transfer mode for the formation of HMXBs. We start by dis- 
cussing the eclipsing ss-HMXB M33 X-7, which hosts a rapidly spinning, 15.65 Mo BH orbiting a close 
to Roche lobe filling underluminous, ~ 70 Mo O-type supergiant in a slightly eccentric orbit with a short 
orbital period of 3.45 days [189]. The existence of this HMXB can neither be explained by CE evolution nor 
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rotational mixing, but dynamically stable non-conservative mass transfer offer a solution [190]. Starting with 
a zero-age MS—MS binary with stellar masses of ~ 90 Mo and ~ 30 Moand an orbital period of ~ 3 days, 
the more massive MS star fills its Roche lobe and dynamically stable mass transfer occurs. A fraction of the 
mass leaving the donor escapes the binary as strong isotropic winds and part is effectively accreted by the 
accretor. After mass transfer finishes, due to the strong winds when the initially more massive star becomes 
a Wolf-Rayet star, the binary hosts a Wolf-Rayet star synchronized with the orbital motion, due to tidal inter- 
action, and a massive O-type MS star. The Wolf-Rayet star with an iron-nickel core later collapses to a BH, 
leading to rapidly spinning BH in a slightly eccentric short-period orbit with a high-mass O-type MS star. 
When the O-type MS star becomes an O-type supergiant almost filling its Roche lobe, the binary becomes a 
BH-HMXB and wind accretion explains its X-ray luminosity. This formation picture is also consistent with 
other sg-HMXBs hosting rapidly rotating BHs and O-type supergiants almost filling their Roche lobes, such 
as Cyg X-1 [191] and LMC X-1 [192], in which dynamically stable non-conservative mass transfer can not 
only explain their masses and orbital periods, but also their BH spins [101]. 

Even though the initial orbital periods in such models are very short (only a few days), they are neither 
odd nor rare, when compared with observations. The orbital period distribution of binaries hosting O-type 
MS stars, the so-called Sana distribution [193], although extending until ~ 10° days, is dominated by systems 
with orbital periods shorter than ~ 10 days. According to the described formation channel, the orbital period 
does not significantly increase after the onset of mass transfer, which explain why sg-HMXBs can be formed 
through dynamically stable non-conservative mass transfer and still have current orbital periods shorter than 
~ 10 days. 

In most known HMXBs, a NS accretes from the circumstellar decretion disk of a rapidly rotating Be star 
[e.g. 50], and their orbital periods are predominantly in range from ~ 10 to ~ 400 days [e.g. 194]. According 
to the wind compression disk model [195], an equatorial disk is formed around a rapidly rotating star due 
to ram pressure confinement by the stellar wind. The disk is formed because the supersonic wind leaving 
the star surface at high latitudes travels along paths that move it down to the equatorial plane, where the 
material passes through a standing oblique shock on the top of the disk. The ram pressure of the polar wind 
thus confines and compress the disk. The emission lines typically observed in an active Be star are attributed 
to this equatorial disk [196]. As mass from this outflowing disk is transferred to the compact object via the 
accretion gate mechanism [52], preferentially during periastron passages or enhanced activity of the Be star, 
the system becomes bright in X-rays, and observable as a Be-HMXB. 

The typical age of a Be-HMXB is ~ 40- 50 Myr [e.g. 197, 198], which is believed to be a direct conse- 
quence of the NS formation and B star mass loss activity time-scales. Given their properties, specially their 
relatively long orbital periods and their rapidly rotating Be stars, they have most likely formed through an 
episode of stable non-conservative mass transfer [e.g. 199, 100], during which the B star spins up. Despite 
the fact that there is only one known Be-HMXB in which the compact object is a BH (MWC 656 [47]), it 
seems more likely that Be-HMXBs hosting BHs are formed during CE evolution [200], although given the 
uncertainties in the stability criterion [201, 125], dynamically stable non-conservative mass transfer is not 
necessarily ruled out. The main problem with having Be-HMXBs hosting BHs formed through CE evolution 
is that the existence of the Be star cannot be easily explained, since a mechanism different from Roche lobe 
overflow would be needed to spin up the B star. Gravitationally focused wind accretion prior CE evolution 
seems to be a good candidate for that [200, 202]. 

In the Corbet’s diagram [203], i.e. the orbital period versus NS spin period diagram, Be-HMXBs exhibit 
a clear correlation, which is explained in terms of the equilibrium period, defined as the period at which the 
rotation velocity at the magnetospheric radius equals the Keplerian velocity [e.g. 204]. If the NS spin period 
is shorter than the equilibrium period, then matter is efficiently ejected from the vicinity of the NS due to the 
propeller mechanism [205], and the NS spin decreases. On the other hand, if the NS spin period is longer than 
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the equilibrium period, then mass and angular momentum is transferred onto the NS thereby increasing its 
rotational velocity. The correlation comes from the fact that the equilibrium period depends on the accretion 
rate, and in turn, on the orbital period. In addition, Be-HMXBs can be separated into two groups in the 
orbital period versus eccentricity diagram, one being those in low-eccentricity orbits (< 0.2) and the other 
one composed of high-eccentricity orbits (= 0.2) [51, 206]. This could then indicate two distinct formation 
channels, with (core-collapse supernova) and without (electron-capture supernova) strong natal kicks during 
the helium star collapse into the NS. 


5.3 Combination of dynamically stable non-conservative mass transfer and 
common-envelope evolution 


The orbits of AM CVns and UCXBs are so close that two episodes of mass transfer must occur to explain 
their orbital periods, at least one of them being CE evolution. Similarly, WR-HMXBs, which are compact 
objects accreting part of the winds from massive helium stars, are also expected to be formed through two 
episodes of mass transfer, given their short orbital periods. However, unlike the previous case, these sys- 
tems can be formed in two episodes of dynamically stable non-conservative mass transfer. There are then 
three possible combinations for AM CVn and UCXBs: (i) double CE evolution, (ii) CE evolution followed 
by dynamically stable non-conservative mass transfer, and (iii) dynamically stable non-conservative mass 
transfer followed by CE evolution; and there is another possible combination for WR-HMXBs, in addition 
to stable non-conservative mass transfer followed by CE evolution, which is (iv) double dynamically stable 
non-conservative mass transfer. In what follows, we will discuss in more detail these formation channels, 
based on results from binary models. 


5.3.1 Evolution through two episodes of common-envelope evolution 


Starting with a zero-age MS-MS binary, if the initial orbital period is not too short nor too long, the 
initially more massive star will fill its Roche lobe after evolving beyond the MS. If the initially more massive 
star is much more massive than its companion, then mass transfer will be dynamically unstable and CE 
evolution commences. The result of such an interaction will be a helium star or a WD, depending on the 
mass and core mass of the donor at the onset of CE evolution, orbiting the initially less massive MS star with 
an orbital period orders of magnitude shorter than initially. In case of a helium star, depending on its mass 
and composition, it will quickly collapse into a WD, NS or BH, in the latter cases this collapse will produce 
a core-collapse or electron-capture supernova. 

If the orbital period of the post-CE binary is not too short, the initially less massive MS will have room 
to grow and evolve off the MS before filling its Roche lobe. If the compact object is much less massive than 
its companion, then another CE will be triggered and the orbital period is shortened even further. In case the 
companion of the compact object is initially a low-/intermediate-mass MS star and if it fills its Roche lobe 
on the FGB, then it will become a helium WD after CE evolution, if its mass is initially smaller than ~ 2 Mo, 
or a helium star if its mass is initially larger than that. On the other hand, if the low or intermediate mass star 
fills its Roche lobe on the AGB, then the compact object companion will likely become a carbon-oxygen 
WD. After the second CE evolution, the post-CE binary will be brought into contact due to emission of 
gravitational waves, becoming in turn an AM CVn or an UCXB, depending on the nature of the compact 
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object, accreting from either a helium/carbon-oxygen WD or a helium star, or even a hybrid WD, in case the 
helium star becomes a WD before filling its Roche lobe. 

This is a standard formation channel for AM CVns and UCXBs [e.g. 102] in which the accreting compact 
object is formed first during the binary evolution. However, even if this scenario ends up having a non- 
negligible contribution to these populations, it cannot account for the existence of some particular types 
of binaries related to AMCVns and UCXBs, such as close detached double helium WDs in which the 
younger WD is more massive than the older WD [e.g. 142], or millisecond pulsars in detached binaries with 
helium WDs [e.g. 207]. Therefore, alternative scenarios involving dynamically stable non-conservative mass 
transfer are needed. 


5.3.2 Dynamically stable non-conservative mass transfer followed by common-envelope evolution 


Starting from a zero-age MS—MS binary, if the orbital period is sufficiently short, the initially more 
massive star will fill its Roche lobe on the Hertzsprung gap, or when it is an unevolved giant. If the initial 
masses of the stars are comparable, or at least the initially more massive star is not much more massive than 
its companion, then mass transfer will be dynamically stable. The orbital period in this case may increase, 
as illustrated in Figure 6, leading then to a binary with an orbital period longer than the orbital period the 
zero-age binary had initially, hosting a (proto) WD and a companion that is more massive than the zero-age 
MS stars. This companion can still be on the MS or already slightly evolved, since it could have started 
its evolution off the MS, in case the initial masses are comparable. The companion of the (proto) WD then 
keeps evolving until it fills its Roche lobe, leading to dynamically unstable mass transfer, since it is much 
more massive than the WD. The post-CE binary will have an orbital period orders of magnitude shorter than 
at the onset of CE evolution, hosting a WD paired with a WD or helium star, which will later collapse into a 
WD, NS or BH, depending on its mass. If the post-CE binary orbital period is sufficiently short, the (proto) 
WD will fill its Roche lobe, giving rise to an AM CVn or UCXBs, depending on the nature of the accreting 
compact object. 

Unlike the double CE evolution scenario, in the scenario in which dynamically stable non-conservative 
mass transfer is followed by CE evolution, the accreting compact object is formed second. This scenario has 
been invoked to explain the existence of double helium WDs in which the younger WD is more massive 
than the older WD [60], which is a natural consequence of this scenario. Other types of double WD binaries 
can also be formed in such a way, e.g. SDSS J1257+5428 [e.g. 208], which is a close binary with orbital 
period of 4.6 hours, hosting a much older extremely-low-mass WD (< 0.3 Mo) orbiting a much more massive 
(~ 1.05 Mo) WD. Most importantly, this scenario can solve the energy budget problem we mentioned earlier 
[e.g. 132, 105], suggesting that the importance of dynamically stable non-conservative mass transfer for 
binary evolution has been underestimated in several contexts. However, the existence of millisecond pulsars 
in detached binaries with helium WDs cannot be explained within this formation channel, which means that 
another channel most likely plays a key role in formation of UCXBs. 

Regarding WR-HMXBs, they are believed to be the direct progenitors of merging binary BHs, which 
means that these merger events can also be used to further constrain WR-HMXBs [e.g. 209]. We have 
already discussed that the properties of ss-HMXBs hosting BHs can be reasonably well explained by an 
episode of dynamically stable non-conservative mass transfer. This is also the most common outcome for 
the first episode of mass transfer in merging BH binary models [e.g. 210, 211]. In these binaries, the BH is 
accreting from the winds of an O-type supergiant that is almost filling its Roche lobe. When it effectively 
fills its Roche lobe, depending on the binary properties, mass transfer can be dynamically unstable or not. A 
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CE evolution is the most common scenario for the formation of merging binary BHs, specially those mergers 
with massive stellar-mass BHs, such as GW190521 [e.g. 212]. 

The sg-HMXB M33 X-7 will apparently undergo CE evolution, given that its O-type supergiant is much 
more massive than its BH companion [189]. On the other hand, in the sg-HMXB Cyg X-1, the O-type 
supergiant is more massive than its BH companion only by a factor of two [191], which implies that mass 
transfer will be most likely dynamically stable when the O-type supergiant fills its Roche lobe. Finally, 
the mass of the O-type supergiant in sg-HMXB LMC X-I is around three times its BH companion mass 
[192], and it is much less clear whether mass transfer could be dynamically stable or unstable in this case. 
Now regarding the current population of confirmed WR-HMXBs, there are two systems that have very 
short orbital periods, namely Cyg X—3, in the Milky Way (4.8 hours) [213], and CG X-1, in the Circinus 
galaxy (7.8 hours) [214], which mostly implies that the BHs orbit within the Wolf-Rayet photospheres. 
These extremely short orbital periods suggest that there was a strong orbital shrinkage, mostly likely only 
possible during CE evolution. For the remaining confirmed WR-HMXBs, the orbital periods are longer than 
a day and these systems could be explained by a weak spiral-in during dynamically stable non-conservative 
mass transfer [53], as we will discuss later. 


5.3.3 Common-envelope evolution followed by dynamically stable non-conservative mass transfer 


The situation here is initially similar to the double CE evolution scenario, i.e. from a zero-age MS—MS 
binary, the initially more massive star fills its Roche lobe when it evolves off the MS, leading to dynamically 
unstable mass transfer, because it is much more massive than its companion. However, unlike in the double 
CE evolution scenario, the post-CE binary orbital period is shorter and the compact object companion is 
not allowed to significantly grow before filling its Roche lobe. Therefore, in this case, the zero-age MS—MS 
binary orbital period must be shorter than that in double CE scenario. 

When the companion of the compact object fills its Roche lobe, depending on its mass and the mass of 
the compact object, mass transfer can be either dynamically unstable or dynamically stable. If the donor is 
much more massive than the compact object, mass transfer will be dynamically unstable, which implies that 
the compact object and its companion will most likely merge on the dynamical time-scale of the donor. In 
case the donor is an MS star, this is the natural outcome of dynamically unstable mass transfer. On the other 
hand, if the donor is a Hertzsprung gap star, or an unevolved FGB star, then a merger is also expected, since 
there would not be enough orbital energy available to eject the CE, due to the short orbital period of the 
binary. 

If the donor is a low-mass star, then mass transfer will most likely be dynamically stable and non- 
conservative. The beginning of the evolution can be thermally stable or not, depending on the masses of 
the accreting compact object and the donor, and on the properties of the donor. In any event, the binary will 
become a CV or a LMXB, depending on the nature of the compact object. If the donor is an M-/K-type MS 
star, i.e. an MS star with mass < 0.8 Mo, then it will eventually become a hydrogen-rich degenerate star, i.e. 
a brown dwarf, when enough mass is lost. However, if the donor is an evolved MS star of an earlier type or 
a Hertzsprung gap star, then it has already or will develop a growing degenerate helium core. In these cases, 
when the donor hydrogen envelope is entirely consumed, it becomes a semi-degenerate helium-rich star. At 
this moment, the CV/LMXB becomes an AM CVn/UCXB. 

The donor transition from a non-generate to a semi-degenerate configuration may be preceded by a short 
detached phase, depending on the post-CE binary properties [e.g. 215, 216]. Alternatively, if the post-CE 
binary properties are such that this detached phase happens at longer orbital periods (= 1 day), then emis- 
sion of gravitational waves will not be strong enough to bring the binary close enough so that the compact 
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object companion could fill again its Roche lobe within a Hubble time. The existence of millisecond pul- 
sars in close detached orbits is essentially explained in such a way, since the NS becomes a millisecond 
pulsar during the LMXB evolution (the so-called recycling scenario [e.g. 217, 218]), and the close detached 
population are those systems observed in the transition from LMXB to UCXB. Those systems that directly 
become UCXBs, or effectively manage to become UCXBs after the detached phase, are then UCXBs hosting 
accreting millisecond pulsars. 


5.3.4 Evolution through two episodes of dynamically stable non-conservative mass transfer 


We have argued before that WR-HMXBs with extremely short orbital periods can most likely only be 
explained with CE evolution. However, the remaining confirmed WR-HMXBs have orbital periods longer 
than ~ | day, namely X-1, in the IC 10 galaxy (34.9 hours) [219], X—1 in the NGC 300 galaxy (32.8 hours) 
[220], and ULX-1, in the galaxy M101 (196.8 hours) [221]. Despite their short orbital periods, they could 
have formed through dynamically stable non-conservative mass transfer, since it is possible in such cases to 
avoid CE evolution during the formation of the Wolf-Rayet stars, due to the high masses of the Wolf-Rayet 
progenitors. 

A slow spiral-in during dynamically stable non-conservative mass transfer is possible if the Wolf-Rayet 
progenitor has a radiative envelope and if it is at most ~ 3 — 4 times the BH mass [222, 223], although the 
exact critical mass ratio depends on the several model assumptions [e.g. 201, 125]. This seems to be the 
fate of several WR-HMXB progenitors. Based on an analysis of the small sample of Wolf-Rayet + O-type 
MS binaries, with orbital periods clustering at around one week [224], it has been recently suggested that, 
after the Wolf-Rayet star becomes a BH and its O-type companion fills its Roche lobe, they are expected to 
become WR-HMXBs with orbital periods as short as ~ 1 day [222]. These WR-HMXBs will then evolve to 
binary BHs, when the Wolf-Rayet donors collapse to BHs. This then suggest that not only the formation of 
WR-HMXBs can avoid CE, but also the formation of merging BHs [201, 125]. 


5.3.5 Further considerations on the formation of ultra-compact X-ray binaries 


The first three scenarios highlighted above correspond to the most traditional scenarios for the formation 
of AM CVn and UCXBs. Regarding UCXBs, there is an additional channel, via accretion-induced collapse 
of an oxygen-neon-magnesium WD [e.g. 225]. If such a WD triggers electron-capture reactions at its centre 
upon reaching the Chandrasekhar mass limit, then the burning that propagates outward makes central tem- 
peratures and pressures high enough to cause a collapse of the WD into a NS [e.g. 226]. Binary models, in 
which different types of donors have been considered, mainly MS stars and helium star donors, have shown 
that this is indeed possible [e.g. 227], leading to the formation of millisecond pulsar binaries. 

In the case of MS stars, after the accretion-induced collapse, the binary temporarily detaches (for 10° — 
10° years), due to the small kicks, but, due to orbital AML, it soon becomes an LMXB. And the subsequent 
formation of the UCXB is similar to that already discussed. The UCXB progenitors in this case are SSXBs 
with MS donors with masses of ~ 2.0—2.6 Mo, for a metallicity of 0.02, and of ~ 1.4 -2.2 Mo, for a 
metallicity of 0.001. It is worth mentioning that this formation channel may also help to explain the existence 
of millisecond pulsars in eccentric detached binaries [e.g. 228]. In the case of helium stars, the binary also 
temporarily detaches and soon becomes an UCXB, due to orbital AML through emission of gravitational 
waves. Unlike the previous case, the helium stars quickly (< 1 Myr) become carbon-oxygen WDs with 
masses of ~ 0.6—0.9, causing the binary to detach. Further orbital AML could eventually make the binary 
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semi-detached again, and the UCXB would host a millisecond pulsar accreting from a carbon-oxygen WD 
donor. This happens for AM CVns with helium star donors with masses of ~ 1.1 — 1.5 Mo 

Even though possible, the contribution from this channel is most likely small. This is because of the nar- 
row range in the parameter space of the accreting WD that allows it to eventually exceed the Chandrasekhar 
mass limit. In particular, during the accreting WD evolution, there is a very narrow range of mass transfer 
rates such that the WD mass can effectively increase [e.g. 12, 229, 230]. In addition, given that a typical 
mass of an oxygen-neon-magnesium WD is ~ 1.2 Mo, the donor cannot be much more massive than that, 
since otherwise mass transfer becomes dynamically unstable, which restricts even more the parameter space 
of the UCXB progenitors. 


5.4 Additional channels through dynamical interactions in high-density environments 


In addition to the formation channels highlighted above, dynamical interactions in crowded regions of star 
clusters, specially globular clusters, can create conditions to form accreting compact objects from systems 
that would otherwise not evolve into them. A way of addressing whether a type of X-ray sources is pre- 
dominantly dynamically formed in globular clusters is by means of the cluster stellar encounter rate, which 
depends on the central density, the core radius and the central velocity dispersion [231, 232]. If dynamical 
interactions are supposed to play any significant role in the formation of X-ray sources, a correlation between 
the number of sources of a particular type within a cluster and its stellar encounter rate could indicate that 
dynamical interactions played a significant role in their formation, despite the several potential problems 
with this interpretation [233]. For a long time, this has been believed to happen for LMXBs and UCXBs, 
which correspond to the brightest X-ray sources found in Milky Way globular clusters, as there is continu- 
ously growing strong evidence for such a correlation [234, 235]. There are currently 21 bright X-ray sources 
in 15 Milky Way globular clusters, being 8 persistent and 13 transients [40]), and they are still consistent 
with being dynamically formed, especially considering their mass density, which is much higher than that 
of systems in the Milky Way disk. 

Given the problem with reproducing the properties of BH-LMXBs in binary models [123] and the above- 
mentioned overabundance of LMXBs in the Milky Way globular cluster population in comparison with the 
Milky Way disk population, a natural solution to the problem could be a dynamical origin for the field popu- 
lation (i.e. those LMXBs not belonging to stellar aggregates). Recent Monte Carlo numerical simulations of 
evolving globular clusters [236] show that progenitors of LMXBs are ejected from their host clusters over a 
wide range of ejection times, contributing to ~ 300 NS-LMXBs and ~ 180 BH-LMXBs to the Milky Way 
field population, not all of them being indeed observable due to observational selection effects. More than 20 
galaxies within ~ 25 Mpc have been deeply investigated with Chandra and Hubble Space Telescope [237]. 
These observations have revealed that, in order to properly describe the X-ray luminosity function of the 
field population of LMXB, a component that scales with the globular cluster specific frequency is needed. 
In addition, it has a shape consistent with that found for the globular cluster population of LMXBs. This 
then indicates that most, if not all, LMXBs in the fields of these galaxies are formed similarly to those in 
globular clusters, and in turn could have origin in globular clusters. These binaries could be ejected from 
the globular clusters due to energetic dynamical interaction leading to their formation, or the cluster itself 
could dissolve due to strong tidal interaction. However, for galaxies with low globular cluster specific fre- 
quencies, the X-ray luminosity function has a component that scales with stellar mass, suggesting then that 
a substantial population of LMXBs in these galaxies are formed through binary evolution, without invoking 
dynamics. 
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These findings may alleviate quite a bit the problems of reproducing the proprieties of BH-LMXBs in the 
Milky Way field. However, it is not clear whether enough systems are predicted to form through this channel, 
nor whether their positions in the Milky Way as well as their kinematic properties, nor their numbers could be 
explained as ejected from globular clusters, or from dissolving star clusters. For instance, while the observed 
systems are almost entirely concentrated in the Milky Way disk [238], those predicted to form in globular 
clusters and latter ejected to the field are spherically distributed around the Milky Way disk, reflecting the 
observed distribution of globular clusters in the Milky Way. In addition, the number of systems predicted 
by this channel seems too low to account for the large number of BH-LMXBs inferred from observations 
(~ 10* — 108 [239]). Therefore, the current theoretical and observational evidence indicates that a dynamical 
origin in star cluster is unlikely to be the dominant formation channel for these systems. A more promising 
alternative to bring into agreement predicted and observed properties seems to be the standard CE evolution, 
coupled with a revision of the strength of orbital AML during their evolution, as we will see later. 

Regarding the fainter X-ray sources such as CVs, the situation is much more difficult to be addressed, due 
to the several biases and selection effects [233], e.g. the photometric incompleteness in the more crowded 
central parts of the clusters is most likely playing a key role in the characterization of these sources [e.g. 240]. 
Therefore, any conclusions drawn from a comparison between the numerical simulations and observations 
of faint X-ray sources should be taken with a grain of salt. Despite that, it is still possible to infer the 
impact of dynamics on this population, if observational selection effects are somehow addressed. Recent 
numerical simulations suggest that fewer CVs should be expected in dense globular clusters relative to the 
Milky Way field, due to the fact that dynamical destruction of CV progenitors is more important in globular 
clusters than dynamical formation of CVs [241]. This is consistent with observations [242], because the faint 
X-ray populations, primarily composed of CVs and chromospherically active binaries, are under-abundant 
in globular clusters with respect to the Solar neighbourhood and Local Group dwarf elliptical galaxies. 
Regarding the dynamical formation of CVs in globular clusters, these simulations show that the detectable 
CV population is predominantly composed of CVs formed via standard CE evolution (= 70%). This happens 
because the main dynamical scenarios have a very low probability of occurring, which resulted in (or very 
weak, if at all) correlation between the predicted number of detectable CVs and the predicted globular 
cluster stellar encounter rate. Despite the fact that these results are consistent with recent observations [242], 
we shall mention that, even though it is not expected to be dominant, dynamics should play a sufficiently 
important role to explain the relatively larger fraction of bright CVs (i.e. those close to the MS in the colour- 
magnitude diagram) observed in denser clusters. 

Regarding SySts and SyXBs, despite a few reasonable candidates have been proposed over the years, 
none so far has been confirmed in the entire Milky Way globular cluster population. This issue has been 
partially addressed in Monte Carlo numerical simulations of evolving globular clusters [164], in which SySts 
with WDs formed without Roche lobe overflow have been investigated. These simulations show that their 
orbital periods are typically sufficiently long (> 10° days) so that, dynamical disruption of their progenitors 
is virtually unavoidable in very dense clusters. On the other hand, in less dense clusters, some of these SySts 
are still predicted to exist, although their identification should be rather difficult in observational campaigns, 
mainly due to their locations within the clusters and rareness. Despite these simulations aimed at explaining 
the lack of SySts in globular clusters, they do not necessarily correspond to definite answers to the issue, 
since most known SySt have orbital period between ~ 200 and ~ 10° days, which means that most SySts 
were not addressed in the simulations. Further simulations are still missing, in which the formation of SySts 
through dynamically stable non-conservative mass transfer is consistently taken into account and the impact 
of dynamics on their progenitors is evaluated. 

So far, we have discussed the population of accreting compact objects in old star clusters. These star 
clusters cannot harbour HMXBs, since at their present-day ages, the entire population of HMXBs have 
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already become single NSs/BHs or binaries hosting either NSs and/or BHs. However, in very young starburst 
galaxies (only a few Myr old), HMXBs are expected to be the dominant source of X-rays. Given that high- 
mass stars, and potentially all stars [243, 244], are most likely formed in binaries [e.g. 193] in embedded 
clusters [e.g. 245], it is not surprising that observations suggest that the HMXBs in starburst galaxies are 
consistent with being formed in star clusters [e.g. 246]. Even though young star clusters are most likely 
HMXB factories, these HMXBs are not necessarily formed dynamically. The analysis of the X-ray sources 
in the merging Antennae galaxies and in the dwarf starburst NGC 4449 suggests that high stellar density, and 
in turn dynamical interactions, are not a primary driver of HMXB formation [247], which is consistent with 
star cluster direct N-body simulations [248]. Therefore, despite being rather frequent in starburst galaxies 
and likely formed in compact star clusters, there is evidence supporting that HMXBs form directly from 
primordial high-mass binaries, without any significant influence of dynamics. 


6 Secular Evolution 


The mass transfer in accreting compact objects is dynamically stable and can be separated into three main 
modes: (i) Roche lobe overflow (if Ra ~ Rp), (ii) atmospheric Roche lobe overflow (if 0.9 RRL $ Ra S RRL), 
and (iii) wind accretion (if Ra < 0.9 RgL). In CVs, LMXBs, AM CVn and UCXBs, the donor is filling its 
Roche lobe. Therefore, mass transfer in these systems can be driven by orbital AML or by nuclear expansion 
of the donor, with or without an initial phase of thermally unstable mass transfer. On the other hand, despite 
the fact that in most SySts, SyXBs, sg-HMXBs and WR-HMXBs, the donor is under-filling its Roche lobe, 
mass transfer still occurs, mainly when the extended atmosphere of the donor fills its Roche lobe, or due to 
accretion of a significant part of the winds from the donor. 

Among the observational properties, the donor mass and the mass transfer rate deserve special attention, 
as these two quantities are somewhat tied to the secular evolution of accreting compact objects. We show 
these properties of accreting compact objects in Figures 7 and 8, against their orbital periods, together with 
evolutionary tracks computed by us with the mesa code and other tracks from previous binary models. In 
what follows, we will describe the dynamically stable evolution of accreting compact objects as well as 
the up-to-date observational information that can be used to constrain evolutionary models, which includes 
their orbital periods, their mass transfer rates, and their donor and compact object masses, but without being 
exclusively limited to them. 


6.1 Cataclysmic variables and low-mass X-ray binaries 


As mentioned earlier, CVs and LMXBs are formed through an episode of CE evolution. More specifically, 
after CE evolution, a post-CE binary is formed in a detached binary hosting a compact object and a non- 
degenerate star. Because of orbital AML and/or nuclear evolution of the future donor, the orbital separation 
decreases up to the formation of the accreting compact binary, when the non-degenerate companion starts 
filling its Roche lobe and becomes a donor. The donors in CVs and LMXBs are usually MS stars, subgiants 
or unevolved red giants. 
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Fig. 7 Comparison between observed and predicted properties of accreting compact objects in the plane donor mass (Ma) 
versus orbital period (Porb). The observational data corresponds to BH-LMXBs [123, 249], BH-HMXBs [220, 192, 189, 47, 
213, 219, 221, 191], NS-LMXBs [250, 78], UCXBs [78], eclipsing ss-HMXBs [251, 252, 253, 254, 255, 256, 257, 258, 259], 
eclipsing AM CVns [260, 261, 262], eclipsing CVs [263], S-type SySts (most of which are eclipsing) [21, 176, 159, 178, 162], 
SyXBs [124, 264], and SSXBs [265]. In particular, we labeled the only known Be-HMXB hosting a BH (MWC 656 [47]), a 
WR-HMXgB with an extremely short orbital period (Cyg X—3 [213]), and the two SySts T CrB [176] and RS Oph [159], which 
are among the few SySts in which the red giant is nearly filling its Roche lobe, and in addition, given their high carbon-oxygen 
WD masses, they will most likely evolve towards Type Ia supernovae [27]. The theoretical evolutionary tracks correspond to 
the lines, of which all coloured tracks were computed by us using the mesa code, adopting solar metallicity, assuming the RVJ 
prescription [73] (Equation 19, with ymp = 3), for the MB, and considering tidal interaction [69] and stellar wind mass loss, 
with the Reimers parameter equals to 0.5 for FGB stars and the Blocker parameters equals to 0.02 for AGB stars. In all these 
coloured tracks, we chose several different initial orbital periods and initial masses for the point-mass compact object and its 
zero-age MS star companion. The black tracks correspond to previous binary evolution models, in which the donor is initially 
a helium WD [266], a helium star [267, 268], or a MS star [265, 215]. Given that the evolutionary tracks pass through most 
observational points, we could naively conclude that our current knowledge on how accreting compact objects evolve is not 
bad. However, major problems arise when we take into account the entire populations of different types of accreting compact 
objects, and confront results from population synthesis with observations. These problems are most likely connected with our 
poor understanding of the mechanisms driving the evolution of accreting compact objects. Despite that, there has been important 
progress in the last couple of years, which can definitely guide us towards a solution to these problems. 


6.1.1 Low-mass unevolved M-/K-type main-sequence star donors 


Let us start with the case in which the donor is initially a low-mass unevolved MS star, corresponding 
to post-CE binaries hosting MS stars with types M and K (< 0.8 Mo), which do not significantly evolve 
within the Hubble time. This case is illustrated in Figures 7 and 8 by the solid red track, for which the 
point-mass compact object and zero-age MS star have both the same initial mass of 0.8 Mo, and the initial 
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Fig. 8 Comparison between observed and predicted properties of accreting compact objects in the plane mass transfer rate (Ma) 
versus orbital period (Porb). We only considered in the plot those systems in which the donor is a Roche lobe filling low-mass 
star, i.e. BH-LMXBs [269], NS-LMXBs [78], UCXBs [78], AM CVns [29], and CVs [77]. Some of the coloured evolutionary 
tracks in Figure 7 are shown, which were computed by us with the mesa code. The black lines correspond to evolutionary 
tracks from previous binary models, in which the donor in initially either a WD [266], or a helium star [267, 268], or a MS star 
(270, 215], the latter having an initial phase of thermally unstable mass transfer. We also included in the figure, as grey thick 
lines, the critical mass transfer rates [269], above which the accretion disks are stable and the systems are persistent sources, and 
transient otherwise, considering irradiated accretion disks or not. By comparing the observational data with these critical mass 
transfer rates, we can see that the disk instability model reasonably well explain persistent and transient sources, except perhaps 
in the case of Cyg X—2 (but see [269]) and some AM CVns. There are two representative tracks for CVs, the solid red line 
computed by us and the long-dashed black line [270], which were both computed with the RVJ prescription [73] (Equation 19, 
with yms = 3). We can clearly see that this MB prescription is not adequate to explain the persistent and transient sources with 
orbital periods longer than ~ 3 hours. In addition, CVs hosting fully convective M-type MS stars, at orbital periods shorter 
than ~ 2 hours, exhibit a spread in the mass transfer rates which cannot be explained by only gravitational wave radiation. 
With respect to NS-LMXBBs, it is also clear that the RVJ prescription (dashed red lines) is not able to account for the persistent 
sources, which has been a long-standing problem of LMXB evolution. However, when a MB law [78] similar to the recently 
developed CARB prescription [79] is adopted (dot-dashed black lines [215]), mass transfer rates are great enough so that they 
can be explained. The persistent UCXBs with long orbital periods (2 40 minutes) can apparently only be explained if the 
donor is initially a helium star, which is consistent with observations [e.g. 267]. There are two particular systems, one AM CVn 
(SDSS J1908+3940) and one CV (SDSS J1538+5123), that are very interesting, since they apparently cannot be explained by 
any of the tracks we included in the figure. While it is has been proposed that this CV might be either a newly born system, or 
a post-nova system, it is not clear what might be happening with this AM CVn. 


orbital period is 12 hours. In this case, the MS donor has a radiative core, which implies that MB is the main 
orbital AML mechanism driving the binary evolution. MB causes mass transfer rates high enough so that 
the MS donor is only able to maintain thermal equilibrium if it is slightly bloated, which means that MS 
donors are somewhat oversized in comparison with isolated MS stars [70]. When the MS donor becomes 
fully convective, MB is expected to cease or, at least, become much less efficient (disrupted MB scenario), 
most likely as a consequence of a rise in the secondary star magnetic complexity [271, 272], rather than 


Contents 39 


the originally proposed demise of the dynamo [273, 274]. From this point on, the main driving mechanism 
becomes GR. This leads to a significant drop in the mass transfer rate and a slowdown of the evolution, as 
AML through GR is in general orders of magnitude weaker than MB. Such a drop in the mass transfer rate 
allows the MS donor to re-establish thermal equilibrium at its normal size (i.e. decrease in size), and the 
system becomes a detached binary since the companion of the compact object is no longer filling its Roche 
lobe. This detachment occurs at the upper edge of the so-called orbital period gap (a paucity of observed 
accreting systems in the orbital period range of 2-3 hrs). 

Even though the system is now detached, it keeps losing angular momentum due to GR and continues to 
evolve towards shorter orbital periods. When the orbital period is sufficiently short, the MS donor fills its 
Roche lobe again, and mass transfer restarts, i.e. the system is again a CV evolving now with lower mass 
transfer rates towards shorter orbital periods, which happens at the lower edge of the orbital period gap. 
When the donor mass becomes too low to sustain hydrogen burning, it becomes a hydrogen-rich degenerate 
object, i.e. a brown dwarf. From this point on, the donor expands in response to the mass loss, leading to 
an increase in the orbital period. The so-called orbital period minimum is then a natural consequence of the 
transition from a fully convective MS star to a brown dwarf. 

The edges of the orbital period gap, and consequently its width, as well as the orbital period minimum 
depend on the metallicity [275, 265], the accreting compact object mass and the orbital AML rate [70, 76]. 
For instance, regarding the impact of the metallicity, the higher the central helium fraction, the shorter the 
orbital period minimum, and the orbital period minimum can be as short as ~ 10 — 40 minutes, when the 
central helium fraction is nearly unity [265]. The orbital period minimum is also affected by orbital AML in 
the following way. For systems having fully convective donors at a given orbital period, the stronger/weaker 
the orbital AML, the higher/lower the mass transfer rate, and in turn the smaller/greater the donor mass, 
which allows the donor to reach the mass limit for hydrogen burning at longer/shorter orbital periods. Or- 
bital AML is also expected to have a strong impact on the orbital period gap, due to variations in the MB 
strength for systems with donors having convective cores. The stronger/weaker MB, the more/less bloated 
the donor is, which means that it becomes fully convective at a lower/higher mass (see Figure 2), having in 
turn a smaller/greater radius. This implies that system will become semi-detached again at a shorter/longer 
orbital period, i.e. the lower edge of the orbital period gap is located at shorter/longer orbital periods. For 
a sufficiently weak MB , such as the RVJ prescription with yms 2 5, the donor bloating is negligible, espe- 
cially near fully convective boundary, and the system evolves virtually without detaching. On the other hand, 
for a very strong MB, such as that provided by the CARB prescription or by the Kawaler prescription with 
n z 1.5, the donor becomes fully convective at a mass as small as ~ 0.18 Mo, and the orbital period gap can 
start at an orbital period as long as ~ 7 — 8 hours. This happens because, as the donor approaches the fully 
convective boundary, it is increasingly driven out of thermal equilibrium, leading to an increase in its radius 
by a factor of = 2 in comparison with isolated M-type MS stars with similar masses. Then, as a response of 
the donor to mass loss, the orbit expands. 

The above outlined evolution successfully describes several features of the evolution of LMXBs and CVs. 
However, it is clearly not complete. For instance, it does not include the impact of the magnetic field of the 
compact object on the evolutionary picture for CV and LMXB, which can be significant as we will explain 
in what follows. If the compact object magnetic moment, which is proportional to its magnetic field strength 
and to the cube of its radius, is sufficiently large (z 10°? G cm’), then a coupling of the donor magnetic 
field lines and the compact object magnetic field lines may lead to an increase of the portion of the donor 
winds that remains trapped to the binary [276, 277]. This effect would make winds to carry away less angular 
momentum, resulting in reduced orbital AML via MB and, in turn, diminished mass transfer rates. While 
NS magnetic fields seem to considerably decay during LMXB evolution [278], the magnetic field strength 
in most accreting WDs are expected to change only negligibly, which is a consequence of mass loss due 
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to nova eruptions [279]. Given the small sizes of NSs (only ~ 10 — 20 kilometers) and their relatively weak 
magnetic fields (~ 10° — 10° G) due to accretion-induced field decay, their magnetic moments are usually 
not strong enough to cause an important impact on LMXB evolution. 

The situation is different for WDs, as they are much bigger (~ 5000 — 10000 kilometers), which implies 
that, if they host sufficiently strong magnetic fields (= 10’ G), the CV evolution is most likely different from 
the remaining CVs hosting WDs with weaker magnetic fields [76]. In case the WD magnetic moment is 
higher than ~ 10°+ — 10° G cm?, MB is completely suppressed and the CV evolution is driven by only GR, 
i.e. at a much lower orbital AML rate, without detaching. Recent binary models provide strong support to 
this magnetic evolutionary picture for such CVs, as the predicted mass transfer rates are in agreement with 
the observations and the orbital period gap is absent in this population, which reasonably well agrees with 
observations [76]. It even seems that CV evolution in general cannot be understood ignoring the origin of 
WD magnetic fields and its impact on orbital AML and exchange between the stellar components and the 
orbit. Including magnetic field generation due to a rotation- and crystallization-driven dynamo [280, 281] 
and the impact of the generated field on the evolution of CVs, can reproduce several observational facts 
which remain otherwise inexplicable, such as the existence of AR Sco, which hosts a radio-pulsing WD 
[282], the relative numbers of magnetic and non-magnetic CVs [1], and the absence of strong magnetic 
WDs among young detached post-CE binaries [283]. 


6.1.2 Subgiant or A-/F-/G-type main-sequence star donors 


Let us now consider the case in which the donor is a nuclear evolved MS star (i.e. near the end of the 
MS) or a subgiant at the onset of mass transfer to the compact object. In other words, we are now looking at 
post-CE binaries hosting MS stars of spectral type G or earlier (= 1 Mo). If the orbital period of the post-CE 
binary is sufficiently long, the future donor star will have enough time to evolve on the MS, or even become 
a subgiant, before filling its Roche lobe. For stars having convective envelopes, the onset of mass transfer 
occurs due to a combination of orbital AML due to MB and nuclear evolution of the future donor, while for 
stars having radiative envelopes, it is triggered only by nuclear evolution of the future donor. The evolution 
for these two types of donors is illustrated in Figures 7 and 8 by the long-dashed and dot-dashed black 
tracks, which were taken from different binary models [270, 265, 215], by the two dashed red and green 
tracks with shortest initial orbital periods (1 day), and by the two yellow tracks (initial orbital periods set 
as | and 3 days). The difference between the three sets of coloured tracks is the initial mass of the compact 
object, which is 0.8 Mo in the red tracks, 1.2 Mo in the green tracks, 1.4 Mo in the yellow tracks; and the 
zero-age MS mass, which is 1.2 Mo for the green and red tracks, and 3.0 Mo in the yellow tracks. 

Depending on the mass and structure of the donor as well as on the mass of the compact object mass, 
the mass transfer can be thermally unstable at the very beginning of the evolution. Thermally unstable mass 
transfer typically occurs when the nuclear evolved MS or subgiant donor is significantly more massive than 
the compact object at the onset of Roche lobe overflow, as in the case of the yellow tracks in Figures 7 and 
8. More generally, for a given donor evolutionary state, the higher the mass ratio at the onset of Roche lobe 
overflow, the greater the chances that the mass transfer will be thermally unstable. The mass transfer rates 
during this phase are typically larger than those provided by thermally stable mass transfer and can exceed 
~ 1077—1076 Mg yr7!, which is high enough to completely drive the donor out of thermal equilibrium. In 
case the compact object is a WD, systems observed at this phase appear typically as SSXBs. The orbital 
period during this thermally unstable phase decreases as a consequence of the donor response to thermal 
time-scale mass loss. However, when mass transfer stabilizes, the donors re-establishes thermal equilibrium 
and the orbital period can further decrease or reverse and increase. 
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Let us now discuss the evolution following the onset of thermally stable mass transfer. For stars with 
radiative envelopes, orbital AML due to MB does not take place, and the orbital period should increase as a 
response to mass transfer. For stars with convective envelopes, there are two possibilities. In case the nuclear 
time-scale of the donor is shorter than its MB-driven mass-loss time-scale, then the binary evolution will be 
driven by the donor nuclear evolution. Mass transfer then proceeds on the nuclear time-scale of the donor, 
the orbital period increases, and the evolution is usually called divergent. On the other hand, in case MB is 
sufficiently strong, such that the MB-driven mass-loss time-scale becomes shorter than nuclear time-scale, 
the evolution will be driven by orbital AML and the orbital period will decrease. In this case, the evolution 
is called convergent and the system will evolve to either a CV or LMXB with a hydrogen-rich degenerate 
donor (in case the MS donor did not have time to start developing a degenerate helium core) or an AM CVn 
or UCXB with helium-rich degenerate donors (in case the donor is initially a subgiant or in case the MS 
donor managed to develop a degenerate helium core). 

The evolution towards a hydrogen-rich degenerate donor is similar to the case in which the donor is an 
unevolved low-mass MS star (compare long-dashed black and solid red tracks in Figure 8). However, the 
evolution towards AM CVns and UCXBs is rather different from that of the previous cases, which can be 
seen by comparing the long-dashed and dot-dashed black tracks with the solid red track in Figure 7, and 
the dot-dashed black track with the long-dashed black and solid red tracks in Figure 8. For instance, the 
orbital period minimum becomes much shorter for donors that were more evolved at the onset of Roche lobe 
overflow. This happens because the donors have helium cores that are partially degenerate, with masses that 
are too low to trigger helium burning. This allows their radii to considerably shrink as the binary evolves, 
so that the binary can reach orbital periods as short as ~ 6 minutes, before the donor become degenerate and 
the orbital period subsequently increases as a response to mass transfer. 

More generally, whether convergent or divergent evolution takes place depends on the strength of orbital 
AML due to MB, on the fraction of mass loss from the binary, on the metallicity and on the binary properties 
at onset of Roche lobe overflow, most importantly the orbital period, donor structure, and compact object 
and donor masses [e.g. 284, 285]. For instance, the impact of the compact object mass can be seen Figures 7 
and 8 by comparing the dashed red and green tracks with the shortest initial orbital period, in which the RVJ 
prescription with yms = 3 was adopted. We can see that the red track in which the point-mass compact object 
mass is 0.8 Mois convergent, while the green track in which the point-mass compact object mass is 1.2 Mois 
divergent. In addition, when the RVJ prescription is compared with the CARB prescription, the bifurcation 
period obtained with the latter is much longer [286, 285, 215, 216]. This is a natural consequence of the 
orbital AML rates provided by the CARB prescription being much larger than those obtained with the RVJ 
prescription. 


6.1.3 Comparison with observations 


Let us start our comparison between theory and observations discussing how good the disk instability 
model [e.g. 4] is to explain persistent and transient systems. By inspecting Figure 8, we can see that in 
general this model successfully explains why some CVs and LMXBs exhibit outbursts, while others do not. 
Irradiation heating of the accretion disk is usually negligible for the stability of disks in CVs [287] while 
in the case of LMXBs X-ray irradiation needs to be taken into account [269]. The net effect of including 
irradiation is that a system at a given orbital period can be persistent for lower mass transfer rates in compar- 
ison with non-irradiated accretion disks which is in general the case (compare CVs and LMXBs in Fig. 8). 
However, the disk instability model apparently fails to explain two systems in Figure 8, namely the CV 
SDSS J1538+5123 and the NS-LMXB Cyg X-2. 
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The time-averaged mass transfer rates in CVs can be estimated with the disk instability model itself 
[287] or by analysing the ultraviolet emission from the WD during quiescence [288], the latter being used 
for SDSS J1538+5123. From the ultraviolet spectra, the WD effective temperature can be estimated, and 
this is considered a good indicator of the mass transfer rate, since the WD is compressionally heated due 
to the accretion process, and the higher the mass transfer rate, the higher the WD effective temperature 
[289, 290]. SDSS J1538+5123 is a dwarf nova (i.e. transient) located below the orbital period gap and 
therefore should according to the disk instability model have mass transfer rates much lower than the ~ 
1.8+0.7 x 107? Mo yr7! inferred from its high WD effective temperature. However, there is no problem 
with the disk instability model here and this high WD effective temperature is most likely related to a recent 
nova eruption [291] and therefore the temperature is not a good indicator of the time-averaged mass transfer 
rate in this particular case. 

The persistent system Cyg X-2 is the only NS-LMXB that does not agree with the predictions of the 
disk instability model, which clearly indicates that this system should be transient. One possibility for this 
disagreement would be that Cyg X-2 is in fact a transient source undergoing a very long outburst [269], 
similarly to what happens with the well-known BH-LMXB GRS 1915+105, which is a transient source and 
has the longest orbital period (33.85 days) of this population. Given that the longer the orbital period, the 
larger the accretion disk, the disks in these wider systems can accumulate enough mass to fuel an outburst 
for decades. The accretion disk in Cyg X-2 is estimated to have a mass of ~ 2.5 x 1078 g, which is massive 
enough to sustain outbursts lasting up to = 80 years[269]. 

Let us move now compare expected secular evolution of CVs and LMXBs with observations. When we 
compare the tracks with the observational data in these figures, we can see that most CVs are consistent 
with having initially unevolved low-mass MS donors (solid red track). However, some CVs clearly reveal 
abnormal (non-solar) abundances in their spectra, indicating that the material has undergone thermonuclear 
processing by the CNO cycle (see [91] for list of systems). These CVs contribute to ~ 5— 15% of the entire 
CV population [292, |], and given that their donors are not expected to possess this material, two scenarios 
have been proposed to explain them. As a first possibility, they are descendants of SSXBs, i.e. they are 
systems hosting initially nuclear evolved MS stars that underwent a short phase of thermally unstable mass 
transfer before reaching their present-day configurations [293, 294]. In this way, the donor has its outer 
layers stripped during this phase, revealing CNO-processed material. The other possibility would be donor 
contamination by the accretion of CNO-processed material from ejected by the WD during nova eruptions 
[295]. 

Even though the first scenario over predicts the number of such CVs by a factor of 2 —6 [293], it has 
been for a long time accepted as most likely. This is because it is apparently rather difficulty for the donor 
to capture enough nova-processed material [91], e.g. because its geometrical cross-section is small and the 
high-velocity material that the secondary geometrically intercepts is not highly nuclear-processed. How- 
ever, it has been recently suggested that these problems can be overcome if a significant part of the ejected 
material leaves the WD with velocities smaller than the escape velocity [91]. This is possible because the 
thermonuclear runaway on the WD causes a strong propagating shock wave, which leaves behind a large 
velocity gradient in the expanding shell. This implies that part of this shell will reach escape velocity but 
the remaining material will be gravitationally bound to the binary, leading to a non-rotating CE-like struc- 
ture surrounding the CV. This nova-induced CE-like structure will then contain a large amount of non-solar 
metallicity material, which can be efficiently accreted by the donor, creating a significant non-solar con- 
vective region. However, given that virtually all CVs are supposed to undergo nova cycles and that only 
~ 10— 15% of them exhibit peculiar abundances, it remains unclear under which conditions this sort of pol- 
lution can efficiently work, and whether it can explain why most CV donors fail in to efficiently accrete. 
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These questions could and should be addressed in future detailed binary models, since their answers are 
potentially related to major problems connected with CV secular evolution, as we will see later. 

Before moving to LMXBs, we just mention a few interesting long period (~ 14—48 hours) CVs that host 
subgiant donors, i.e. BV Cen [296] and GK Per [297] in the Milky Way field, and AKO 9 in the globular 
cluster 47 Tuc [298]. Interestingly, the masses of the WDs in BV Cen (~ 1.02 — 1.44 Mo [299]) and GK Per 
(~ 0.92 — 1.19 [297]) are significantly higher than the average among CVs (~ 0.840.2 Mo [300, 263, 77]), 
which is consistent with the evolutionary picture in which the WD mass has grown due to the high mass 
transfer rates expected when the donor is initially an evolved MS star or a subgiant. This makes this type of 
CVs very good candidates for being Type Ia supernova progenitors, under certain conditions, although they 
correspond to a negligible fraction of the entire CV population. 

When we compare the distribution of LMXBs and CVs in Figure 7, together with the evolutionary tracks, 
it very likely that the contribution of systems with initially evolved donors is much larger for LMXBs in 
comparison to CVs, as they have on average significantly longer orbital periods for the same range of donor 
masses. While the evolutionary sequence given by solid red line in this figure seems to nicely explain at least 
around half of NS-LMXBs, a lot of them can only be explained if the donor is initially an evolved MS star 
or a subgiant. The situation is even more dramatic for BH-LMXBs, since only two systems host M-type MS 
donors, and the majority of systems with orbital periods shorter than ~ 1 day have K-type MS donors [249]. 
In addition, it seems that a significant fraction of LMXBs start their evolution having initially intermediate- 
mass donors [301], especially in the case of BH-LMXBs. Similarly to CVs, the LMXBs with longest orbital 
periods (= 1 days) most likely have subgiant or red giant donors, such as Cyg X—2 [302], GRO J1744—28 
[303], GX 13+1 [304], XTE J1550-564 [305], GX 339-4 [306], GS 1354—64 [307], GRO J1655—40 [308], 
V4641 Sgr [309], GS 2023+338 [310], and GRS 1915+105 [311]. These systems are then evolving towards 
longer orbital periods, as mass transfer is driven by the nuclear expansion of the donor. 

Even though the evolutionary picture that we have discussed reasonably well explains the main observa- 
tional features of individual CVs and LMXBs, as illustrated in Figures 7 and 8, there are still major problems 
when we take into account the entire populations and compare observed properties with those predicted by 
population synthesis which are most likely intrinsically connected with the orbital AML mechanisms driv- 
ing the binary evolution. In particular, the standard RVJ prescription (Equation 19) for MB usually assumed 
for models of CVs and LMXBs is probably the origin of most of the problems we are currently facing. For 
instance, the predicted fraction of CVs hosting brown dwarfs is much higher than observed [76, 1], which 
suggests that the overall evolutionary time-scale should be much longer than predicted by models using 
the RVJ prescription. In addition, the predicted mass transfer rates of CVs above the orbital period gap are 
typically much higher than those inferred from observations, except for systems near the upper edge of the 
orbital period gap, for which predicted rates are much smaller than those inferred from observations [76, 77]. 
This suggests that the dependence of MB with the orbital period should be the opposite of what is predicted 
by the RVJ prescription. While this recipe leads to a decreasing orbital AML rate as the CV evolves (see 
Figure 2), observations suggest that the overall rate should increase as the CV evolves. Moreover, the pre- 
dicted mass transfer rates for NS-LMXB are so low that they cannot explain the persistent sources [74], 
which suggests that MB should be much stronger in these systems than predicted by the RVJ prescription. 
Furthermore, the predicted donor mass distribution of BH-LMXBs peaks at a much higher value than the 
observed distribution [123], which suggests that MB should also be stronger in these system. Finally, when 
the RVJ prescription is adopted, a well-known fine-tuning problem arises when trying to reproduce millisec- 
ond pulsars in compact binaries (orbital periods in the range ~ 2 —9 hours) with helium WD companions of 
mass $ 0.2 Mo, [286, 215], which represents a big problem in the formation of these binaries. 

Regarding LMXBs and their descendants, recent binary models [78, 79], in which the CARB prescription 
has been proposed, indicate that the problem with the persistent NS-LMXBs can be solved with this recipe 
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which provides much higher orbital AML rates (see Figure 2) and therefore higher mass transfer rates, 
especially for persistent systems with shorter orbital periods. For this reason, the CARB prescription might 
also solve the problem with BH-LMXBs, since it could shift the predicted donor mass distribution towards 
lower values and hopefully reproduce the observed distribution. In addition, this new prescription might 
also be consistent with the fast orbital decay observed in some BH-LMXBs (e.g. XTE J1118+480 [312], 
A0620-00 [313], and GRS 1124-68 [75]), which cannot be explained with the RVJ prescription. When 
applied to compact millisecond pulsar binaries, the CARB prescription also provides a solution to the fine- 
tuning problem [216]. In addition, not only compact binaries are well-reproduced by this prescription but 
also wide-orbit binaries can be explained, and in general, the relation between the helium WD mass and the 
orbital period predicted by this prescription is in agreement with observations [216]. 

Regarding CVs, it has been proposed that not only the problem with the mass transfer rates but also the 
problem with the fraction of post-orbital period minimum could be solved with another MB prescription 
[76, 1]. A reasonable candidate to solve both problems is the Kawaler prescription with n ~ 1.25, which is 
shown in Figure 2. This prescription predicts an overall evolutionary time-scale substantially longer than 
that obtained with the RVJ prescription, which may perhaps contribute significantly to reduce the predicted 
number of CVs that manage to reach the orbital period minimum. Moreover, given that the orbital AML rates 
provided by this prescription are around one order of magnitude higher than GR, the constraints provided by 
detached WD-MS post-CE binaries [314] are not expected to be violated, which strengthens even more the 
potential of solving CV problems with this prescription. 

In addition to the problems most likely connected with the standard RVJ prescription for MB, there 
are other problems that are probably associated with another orbital AML mechanism, i.e. consequential 
AML which is not included in the standard evolutionary picture. For instance, the observed CV WD mass 
distribution has a peak at ~ 0.8 Mo [300, 263, 77], while most predicted CV WD masses are on average 
smaller than that [300, 263, 77, 315]. In addition, for CVs hosting fully convective MS donors, there is a 
clear spread in their mass transfer rates inferred from observation which is definitely not predicted when 
GR is taken as the only driver of their evolution [77]. Moreover, still when only GR is considered, the 
predicted mass transfer rates among post-orbital-period-minimum CVs are much lower than those inferred 
from observations [77], and the predicted orbital period minimum is significantly shorter than observed 
[316, 317, 70]. These problems might be solved if the total AML for these systems is stronger than predicted 
by GR [70, 90, 318]. In case there is another orbital AML mechanism driving the evolution together with 
GR for CVs with fully convective MS donors and brown dwarf donors, then the predicted mass transfer rates 
might be significantly enhanced, which would cause the CVs to evolve faster so that their donors become 
degenerate at longer orbital periods. The problem with the WD masses can be solved if the strength of this 
additional orbital AML is inversely proportional to the WD mass, such that most CVs with low-mass WDs 
would become dynamically unstable at the onset of mass transfer [90]. Similarly, the spread in the mass 
transfer rates for CVs with fully convective MS donors might also be solved with this dependence on the 
WD mass, if perhaps the additional orbital AML also depends on the age of the CV. 

The most natural choice of additional orbital AML to enhance the mass transfer rates would be conse- 
quential AML, and there are very good candidates for the above-highlighted dependence on the WD mass, 
and eventually on the time since the onset of mass transfer[89, 90, 318]. The frictional orbital AML produced 
by novae depends strongly on the expansion velocity of the shell [89]. For low-mass WDs (< 0.55 Mo), the 
expansion velocity is small [11], and this may lead to strong orbital AML by friction that makes most CVs 
with low-mass WDs dynamically unstable. In particular, given the velocity gradient in the nova shell we 
mentioned earlier, it is reasonable to expected that a nova-induced CE-like structure surrounding the binary 
will be formed after the eruption. This idea has been tested in binary models [318], which confirms that 
such a configuration drastically affects the stability of mass transfer in CVs with low-mass WDs. As we 
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mentioned earlier, this form of consequential AML could also explain the existence of CVs with non-solar 
abundance, as an alternative to the post-SSXB scenario, which predicts more such CVs (a factor of 2 — 6) 
than observed. Other promising mechanisms are circumbinary disks [86, 93], or a combination of both, e.g. 
if the circumbinary disk is fed by a fraction of the accreted matter that is re-emitted by WD due to the nova 
eruptions. 


6.2 AM CVns and ultra-compact X-ray binaries 


As we mentioned earlier, AM CVns and UCXBs can be formed in three ways, either through two CE evo- 
lution, or through dynamically stable non-conservative mass transfer followed by CE evolution, or through 
CE evolution followed by dynamically stable non-conservative mass transfer. Since the later channel from 
CVs/LMXBs has just been addressed, we focus on the other two, which have always a CE evolution for 
the second episode of mass transfer. Similarly to CVs and LMXBs, after CE evolution, a detached post-CE 
binary is formed and evolves towards shorter orbital periods due to orbital AML, leading to the onset of mass 
transfer, when the donor fills its Roche lobe. However, unlike these systems, the main mechanism driving 
the evolution of AM CVns and UCXBs is always GR, since typically their donors are initially helium WDs 
or helium stars. 


6.2.1 Helium white dwarf or helium star donors 


Let us start with the case in which the donor is initially a helium WD. This case is illustrated in Figures 7 
and 8 by the dotted black track in each figure, which were taken from a binary model [266]. Even though 
the orbital period corresponding to the onset of mass transfer depends on the mass ratio and on helium WD 
properties, mainly its radius, and in turn its mass, it is of order of a few minutes. After that, the further binary 
evolution depends on the initial entropy and the subsequent thermal evolution of the donor, and is usually 
separated into three phases [e.g. 319, 102, 320, 266]. 

In the very beginning, the donor contracts in response to mass loss as the outermost radiative layer is 
stripped off, leading to an increase in the mass transfer rate. The mass transfer rate then grows from zero to 
its maximum, which can be up to ~ 1078 — 1076 Mo yr7!, depending the donor properties, while the donor 
radius decreases to a minimum value. During this phase, the mass loss time-scale becomes much shorter than 
the donor thermal time-scale, and the orbital period decreases as a response to mass transfer. After this quick 
phase, which lasts < 10° years, the donor radius evolution reverses and the donor expands in response to 
mass loss. This happens because, during the first phase, the entropy profile eventually becomes sufficiently 
shallow that the expansion of the underlying layers starts to dominate. The binary then evolves towards 
longer orbital periods driven by GR, with decreasing mass transfer rates. During this second phase, which 
can last ~ 1 Gyr, the mass loss time-scale is much shorter than the donor time-scale and the donor responds 
adiabatically to the mass loss. The third phase starts when the mass loss time-scale becomes comparable to 
or shorter than the donor thermal time-scale in the underlying layers. This allows the donor to cool down 
and contract, although the contraction ends when the donor has shed sufficient entropy to reach its fully 
degenerate configuration. This phase starts at an orbital period of ~ 40 minutes. 

Let us now address the case in which the companion of the compact object has not developed a degenerate 
core before CE evolution, resulting in turn in a post-CE binary hosting the compact object and a helium star. 
In case the post-CE binary orbital period is sufficiently short, the compact object companion will still be 
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a helium star, i.e. still burning helium, at the onset of mass transfer. The binary evolution for this case 
is illustrated in Figures 7 and 8 by the solid and short-dashed black tracks taken from two binary models 
[267, 268]. Depending on the mass, the lifetime of a helium star are typically ~ 100— 10° Myr, more massive 
stars having shorter lifetimes, and its radius can increase up to ~ 30% during its life [321]. Thus, unlike the 
previous case, the onset of Roche lobe filling is driven by a combination of orbital AML driven by GR and 
donor nuclear expansion. Therefore, the orbital period at which the onset of mass transfer occurs depends 
on the masses of the compact object and helium star as well as on the helium star evolutionary status. In 
particular, the longer the post-CE binary orbital period, the more time the helium star has to evolve, becoming 
larger, and in turn the longer the orbital period at the onset of mass transfer. Given the orbital AML rates 
provided by GR and the helium star lifetimes, to have a helium star at onset of mass transfer, the post-CE 
binary orbital must be typically shorter than ~ 2 hours [321]. After mass transfer starts the further evolution 
for systems hosting helium star donors is different from those having helium WDs [322, 321, 102, 268]. 

During the Roche lobe overflow, the mass loss time-scale is typically shorter than the thermal time-scale 
of the donor which makes the donor slightly bloated as a response to mass loss. This causes a rapid increase 
in the mass transfer rate, reaching values of a few 1078 Mo yr7!. Since the donor should be modestly 
evolving and expanding, in the absence of orbital AML the orbital period should increase in response to 
mass loss. However, since GR is sufficiently strong, the donor nuclear evolution time-scale is longer than 
the GR-driven mass loss time-scale, which leads to convergent evolution, and the orbital period decreases. 
A significant fraction of the energy released through nuclear burning is absorbed in the donor envelope, 
which leads to a decrease in its luminosity and an increase of its thermal time-scale. As the binary evolves 
and the donor mass drops, nuclear burning quickly becomes less important, which makes the convective 
core to be replaced by an outer convection zone that penetrates inwards. At the same time, as the orbital 
period decreases, the mass loss time-scale also decreases, causing the mass transfer rate to slightly increase. 
Towards the end of this initial phase (< 1—5 Myr), the donor thermal time-scale becomes much longer than 
the mass loss time-scale, and the thermal structure of the donor is severely perturbed. The mass transfer 
rate then increases faster, as the donor tries to remain within its Roche lobe. At some point, helium burning 
in the centre of the donor virtually vanishes and the donor is driven far out of thermal equilibrium. The 
mass transfer rate then reaches its maximum and the orbital period its minimum. This typically occurs when 
the donor mass is ~ 0.20—0.26 Mo, at an orbital period of ~ 8—11 minutes. Since the star becomes more 
degenerate, its radius starts to increase and the further evolution resembles that of systems hosting helium 
WD donors, in which the donor is evolving adiabatically in response to mass loss. The mass transfer rate 
drops while the binary evolves towards longer orbital periods. 

In case the compact object is a WD, some interesting features are expected in the first phases of AM CVn 
evolution hosting initially helium star donors. It has been shown by binary models [323, 324] that, if the 
initial mass of the helium star is = 0.4 Mo and the accreting WD mass is initially = 0.8 Mo, enough mass is 
accumulated onto the accreting WD, leading to a first thermonuclear flash that is likely vigorous enough to 
trigger a detonation in the helium layer. This thermonuclear runaway originated from the detonation in the 
helium layer of the accreting WD may create converging shock waves that reach the core. If the carbon in 
the core is also detonated, the thermonuclear runaway is seen as a Type Ia supernova. Those that survive the 
first flash and eject mass will have a temporary increase in orbital separation, but orbital AML through GR 
drives the donor back into contact, resuming mass transfer and triggering several subsequent weaker flashes. 

Despite the fact that the three evolutionary channels leading to AM CVns and UCXBs are in principle 
viable, there are theoretical arguments against two of them. It has been argued that most, if not all, WDs 
accreting initially from helium WDs merge at the beginning of the evolution [325]. This would happen due 
to nova eruptions during the initial phases of accretion which would trigger dynamical friction within the 
expanding nova shell, shrinking the orbit and causing the mass transfer rate to become dynamically unstable. 
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If that is true, most, if not all, AM CVns would come from WDs initially accreting from helium stars or from 
CVs with subgiant or nuclear evolved MS donors. It has been argued, based on the RVJ prescription for 
orbital AML due to MB, that only a negligible fraction of AM CVns could have originated from CVs [270, 
326]. The reason for that is the narrow range of CV properties from which AM CVn originate. However, the 
situation here is very similar to the fine-tuning problem related to the formation of close detached millisecond 
pulsar binaries. As the CARB prescription, which predicts higher orbital AML rates than the standard RVJ 
prescription, can solve the millisecond pulsar fine-tune problem, it would not be surprising if it can also 
increase the probability of forming AM CVns from CVs. 


6.2.2 Comparison with observations 


After describing the main features of AM CVn and UCXB evolution, we can confront prediction with ob- 
servations. Let us start with a brief discussion of the persistent and transient sources with respect to what the 
disk instability model predicts. By inspecting Figure 8, we can see that this model provides quite convincing 
results, similarly to the case of CVs and LMXBs, especially for UCXBs [267]. However, some AM CVns 
deserve attention, like ASASSN-14cc and ASASSN-14mv. The disk instability model predicts that these 
two systems should be persistent, but they have undergone outbursts [29]. In addition, SDSS J1525+3600 
and GP Com have never shown outbursts [29], although the model predicts that they should be transient. Fi- 
nally, despite their very low mass transfer rates ($ 1071! Mo yr7!), SDSS J1208+3550, SDSS J1137+4054, 
V396 Hya and SDSS J1319+5915 are marked in the figure as persistent, since they have never been observed 
in outburst. Even though it is very likely that the lack of outbursts related to the last four systems is due to 
the fact that their accretion disks are stable and cold, it is not clear what might be happening with the other 
systems. 

When we compare the evolutionary tracks for the three channels leading to AM CVns and UCXBs in 
Figures 7 and 8, we can see that they can reasonably well explain observations, and some particular systems 
can be better described by particular channels. For instance, as noted in previous studies [267], UCXBs can 
be separated into three groups: (i) transients with orbital periods longer than ~ 40 minutes, (ii) persistent 
with orbital periods shorter than ~ 25 minutes, and (iii) transient with orbital period longer than that. Since 
the mass transfer rate drops as an UCXB evolves with a WD donor, sources in the group (i) and (ii) could 
in principle be explained by any scenario. However, those in the group (iii), i.e. the long-period persistent 
UCXBs 4U 1626-67, 4U 1916-053, and 4U 0614+091, can apparently only be explained if the donor is 
initially a helium star or perhaps if they are descendants of LMXBs. This is because if the donor is initially 
a helium WD, mass transfer rates in this orbital period range are expected to be similar to those in group (i), 
i.e. orders of magnitude lower than those inferred for the group (iii). 

This has important implications when we compare the population of UCXBs in the Milky Way field with 
those belonging to globular clusters. As already noted before [327, 267], the persistent UCXBs in globular 
clusters have orbital periods much shorter than those of 4U 1626-67, 4U 1916-053, and 4U 0614+091 
(the longest orbital period is 22.58 minutes for M15 X—2 [328]). This fact can be naturally explained if 
long-period persistent UCXBs have initially helium star donors [267]. This is because UCXBs in globular 
clusters cannot be formed in such a way, and the reason for that is twofold. First, the typical evolutionary 
time-scale for an UCXB to be seen as persistent is only a few tens of Myr [268]. Second, only MS stars 
with masses 2 2 Mo will evolve without developing a degenerate helium core on the FGB, which is a much 
higher mass than that of MS stars close to the turn-off point in globular clusters. This would then explain the 
lack of long-period persistent UCXBs in globular clusters. 
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There is one AMCVn (SDSS J1908+3940) that seems challenging to be explained by any evolution- 
ary channel. This system has an orbital period of ~ 18 minutes and an inferred mass transfer rate of 
~4x1077 -107% Mo yr7! [29]. For its orbital period, this mass transfer rate is orders of magnitude higher 
than expected from any type of initial donor, which is $ 1078 Mo yr7!. Future observational efforts and 
binary models with focus on this system could help to understand whether any of these formation channels 
could explain its properties. If not, perhaps something atypical is occurring with this system which pushed 
its inferred mass transfer rate to such a high value. 

Unfortunately, except for a few particular systems, we still do not know which of these three formation 
channels is the dominant among AM CVns and UCXBs, if any is dominating at all, which compromises 
studies of their intrinsic populations. It seems clear as well that, for the overwhelming majority of systems, 
having only the orbital period, donor mass and mass transfer rate is not enough to distinguish between the 
different evolutionary scenarios, and more information is needed. There are a few ways to better evaluate the 
contribution of each channel, by gathering more information about these systems. 

For instance, by identifying and studying the properties of potential progenitors, such as the system 
LAMOST J0140355+392651 [329]. This system has an orbital period of 3.81 hours and hosts a bloated, 
relatively cool, low-mass (~ 0.15 Mo) proto WD and a massive (~ 0.95 Mo) WD companion. Its properties 
can be well-explained as a CV entering the detached phase, when the donor envelope has been entirely 
consumed. In this case, the proto WD was a nuclear evolved MS star before the onset of mass transfer. 
The system will evolve through this detached phase due to GR and mass transfer will eventually resume, 
after a few Gyr. This sort of investigation provides key information to understand binary evolution leading 
to AMCVns and UCXBs. However, the number of these systems are far too low to provide conclusive 
evidences for a particular channel being dominant or not. 

Another possibility is by constraining the properties of their donors with detailed analyses of their chemi- 
cal composition, which allows us to investigate in more detail their evolutionary history [103]. Those formed 
through CV/LMXB evolution are expected to retain a greater fraction of their hydrogen in comparison with 
systems formed in the other two channels. This would imply that the presence of hydrogen in their spec- 
tra is needed to rule out the other channels as well as systems not exhibiting hydrogen in their spectra are 
most likely not formed in such a way. In both helium WD and helium star channels, the transferred ma- 
terial is expected to be mainly composed of helium with CNO-processed material, mainly nitrogen, and 
abundance ratios of these elements are believed to be an effective discriminant between these two channels. 
A third possibility is by means of the donor mass-radius measurements obtained from precise analyses of 
eclipsing systems, since each evolutionary channel predict on average different radius for the same mass 
[319, 321, 270], although some overlap is also expected depending on the initial properties of the donor as 
well as how the WD donor cooling is treated. In general, systems having initially helium WD donors evolve 
with smaller radii than systems having initially helium star donors, and these evolve with smaller radii than 
those originated from CVs/LMXBs. 

The last option we discuss is applicable to UCXBs and consists of investigating the existence/lack of 
thermonuclear explosions that occur on the surface of the NS as well as their features. It has been proposed 
that there is a strong correlation between the iron Ka line strength in their X-ray spectra and the abundance of 
carbon-oxygen in the accretion disc [330]. This line occurs at ~ 6.4—6.9 keV and is a result of the reflection 
of the radiation by the accretion disk and the surface of the WD facing the NS, and for this reason it can 
provide clues on the disk chemical composition. This correlation coupled with the fact that UCXBs with 
prominent and persistent iron Ka emission also show bursting activity, it has been proposed that UCXBs 
with persistent iron emission have helium-rich donors, while those that do not, likely have carbon/oxygen 
or oxygen-neon-magnesium-rich donors [42]. This possibility can be further explored, but more elaborated 
X-ray reflection models will be eventually needed. 
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Having in mind these possibilities, we illustrate with eclipsing AM CVns how they may help to disentan- 
gle the formation channels. Thus far, there are seven eclipsing AM CVns with accurate measurements of the 
masses and radii of their donors [260, 261, 262]. When considering only the radii and masses, YZ LMi can be 
easily explained by either the helium WD or the helium star channels, ZTF JO407—0007 by the helium WD 
channel, and ZTF J2252-0519 by the helium star channel. The remaining (Gaial4aae, ZTF J1637+4917, 
ZTF J0003+1404 and ZTF JO220+2141) have very large radii, which could be only explained by the CV 
channel. However, when the chemical properties of these seven systems are taken into account, the situa- 
tion becomes much more complicated. Regarding those systems apparently formed through the CV channel, 
since no hydrogen was detected in the spectra of any of them, it seems more likely that they have formed in 
a different way. Therefore, the situation for these systems apparently seems rather challenging. A possibility 
to bring into agreement predicted and observed WD donor radii and chemical properties is by avoiding the 
phase in which the donor can contract and cool down. This phase is predicted to exist because, at some point, 
the donor thermal time-scale becomes comparable or shorter than the mass loss time-scale, since the orbital 
AML rate is continuously dropping as the AM CVn evolves. If there is an additional source of orbital AML, 
e.g. consequential AML, such that the donor cooling is significantly delayed, then the donors could keep 
their large sizes, even when the orbital period is longer than ~ 40 minutes. However, this possibility at the 
moment is very speculative and future detailed binary models could verify whether consequential AML is 
able to do that, and if so, for which mechanisms and under which conditions. 

After describing compact objects accreting from Roche lobe filling donors, we can turn to a discussion 
of other modes of transfer, in which the donor is not filling the Roche lobe and is either a low-/intermediate 
mass evolved red giant, or a supergiant, or a Wolf-Rayet star. 


6.3 Symbiotic Stars and Symbiotic X-ray binaries 


We have argued before that, in order to explain their orbital periods (~ 200 — 10° days), most S-type SySts 
should form through an episode of dynamically stable non-conservative mass transfer. After this event, the 
resulting post-stable-mass-transfer binary, with an orbital period longer than ~ 100 days, hosts a WD orbiting 
a low-mass MS star or perhaps a subgiant or red giant, depending on the mass ratio of the zero-age MS—MS 
binary. On the other hand, those S-type SySts hosting massive WDs as well as SyXBs are most likely formed 
through CE evolution, leading to long-period post-CE binaries, since dynamically stable non-conservative 
mass transfer cannot explain their properties. Finally, given the very long orbital periods of D-type SySts 
(= 40 years), their WD progenitors are not expected to form during an episode of Roche lobe filling, i.e. their 
WDs are formed similarly to single WDs. In the three situations, given the relatively long orbital period of 
the binary, the WD companion has enough room to evolve and expand, and the binary will become a SySt 
or a SyXB under certain conditions. Even though in what follows we discuss the evolution of these systems 
with focus on SySts, the case of SyXBs is very similar, except from the fact that the compact object is a NS. 

The donors in SySts are either located near the tip of the FGB or on the AGB [e.g. 331, 22]. The typical 
luminosity of accreting WDs in SySts is ~ 10? — 10+ Lo [23], and in most cases these are powered by 
nuclear burning due to the very high accretion rates, which are typically ~ 1078 — 1077 Mo yr™!. In addition, 
in most SySts, the red giant is under-filling its Roche lobe filling [24, 332], but in some cases it is nearly 
filling its Roche lobe, as evidenced by the clear signature of ellipsoidal variations in their light curve [331]. 
Therefore, in what follows, we discuss two promising modes of mass transfer for these systems, which are (i) 
atmospheric Roche lobe overflow [55], if 0.9RrL_ < Ra S Rex, or (ii) gravitationally focused wind accretion 
[170, 171, 172, 173, 174, 175], if Ra $ 0.9 Rev. 
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6.3.1 Atmospheric Roche lobe overflow 


As the red giant evolves, its radius approaches its Roche lobe and at some point the mass transfer rate 
provided by Equation | will become relevant. If the red giant mass is not much higher than the compact 
object mass, Zaq > ¢RL, and mass transfer is dynamically stable. We show in Figure 7, as dashed green lines, 
several evolutionary tracks in which the initial point-mass compact object mass is 1.2 Mo, and the zero-age 
MS mass is 1.25. For those with initial orbital periods longer than ~ 5 days but shorter than ~ 100 days, the 
onset of mass transfer occurs when the donor is an unevolved red giant on the FGB. The mass transfer in this 
case is similar to what is depicted in Figure 6, and it is characterized by an initial phase of thermally unstable 
mass transfer, until the donor mass is sufficiently low such that the mass transfer stabilizes. After that, mass 
transfer proceeds on the nuclear time-scale of the donor. The mass transfer rates nicely correlate with the 
evolutionary status of the donor. First, the longer the initial orbital period, the more evolved the donor at 
the onset of mass transfer. Consequently, the more evolved the donor, the higher the mass transfer rates. For 
instance, during thermally unstable mass transfer, the rates for the binaries with shortest and longest initial 
orbital period are ~ 3 x 1078 and ~ 6 x 1076 Mo yr™!, respectively. During nuclear time-scale mass transfer, 
the rates for the same binaries are ~ 5 x 107° and ~ 7x 1078 Mo yr™!, respectively. 

There are at least three SySts in which their red giant donors are very close to filling their Roche lobes 
and their WDs are massive enough so that mass transfer can be regarded as dynamically stable. T CrB is 
a well-studied SySt with a relatively short orbital period (227 days), hosting a massive carbon-oxygen WD 
(1.37+0.13 Mo) and a red giant with mass 1.12+0.23 Mo [176]. RS Oph is another well-known SySt, which 
also host a massive carbon-oxygen WD (1.2— 1.4 Mo). However, its orbital period (453 days) is longer than 
that of T CrB and it hosts a less massive red giant donor (0.68 — 0.80 Mo) [159]. This likely indicates that 
both SySts could share similar evolutionary sequences, in which the orbital period is expanding as the red 
giant is consumed. The last SySt we know in which the red giant donor is nearly filling its Roche lobe 
is V3890 Sgr [178]. Similarly to the other two systems, this SySt also host a massive carbon-oxygen WD 
(1.35+0.13 Mo). However, its longer orbital period (747 days) combined with the mass of its red giant donor 
(1.05 + 0.11 Mo) suggests that the evolutionary pathways of V3890 Sgr is slightly different from those of 
T CrB and RS Oph, starting at a longer orbital period. 

The properties of T CrB and RS Oph can be nicely reproduced if the post-CE binary hosts a WD with 
mass ~ 1.15 Mo, and a MS star companion with mass ~ 1.4 Mo, as illustrated in Figure 7 by the solid green 
track passing through these two systems. As mentioned earlier, T CrB and RS Oph could indeed be explained 
with the same evolutionary sequence, if the WD mass in the former is smaller, which would imply that T CrB 
may evolve to something similar to RS Oph. In order to reproduce the properties of these two systems with 
the same evolutionary sequence, the post-CE binary orbital period needs to be as long as ~ 200 days, and 
the fraction of mass lost from the vicinity of the WD due to isotropic re-emission needs to be By ~ 0.6. The 
predicted masses of the WD and red giant during the T CrB stage are ~ 1.25 and ~ 1.07 Mo, and during the 
RS Oph stage they are ~ 1.38 and ~ 0.67 Mo. The WD will exceed the Chandrasekhar limit and explode in 
a Type Ia supernova when the orbital period is ~ 500 — 800 days, depending on the assumed Chandrasekhar 
limit. After the onset of mass transfer, the binary reaches the T CrB stage after ~ 0.2 Myr. The transition 
from the T CrB stage to the RS Oph stage takes ~ 1 Myr, and from RS Oph stage to the Type Ia supernova 
< 1-4 Myr, depending on the assumed Chandrasekhar limit. 

The properties of V3890 Sgr can also be reproduced with similar WD and companion masses after CE 
evolution, as well as the same mass loss fraction during SySt evolution. However, the post-CE binary orbital 
period must be much longer (~ 600 days) than in the previous case. An example of evolutionary sequence 
(rightmost solid green track) is shown in Figure 7. The predicted masses of the WD and red giant at the 
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orbital period of V3890 Sgr are ~ 1.25 and ~ 0.96 Mo. Similarly to the previous case, the evolution will 
terminate when the WD exceeds the Chandrasekhar limit. 

As mentioned earlier, mass transfer will be dynamically stable only if the red giant mass is comparable to 
or only slightly more massive than the compact object mass. This condition is usually not satisfied in SySts 
[21], since their WD masses are typically < 0.6 Mo, while their red giant donors usually have masses of 
~ 1.2-2.2 Mo, implying mass ratios 2 2.5. Such mass ratios are likely high enough to cause mass transfer 
to be dynamically unstable. In this case, as the red giant radius approaches the Roche lobe radius, the mass 
transfer rates provided by Equation | always increase, leading to CE evolution. This behaviour is illustrated 
in Figure 9 by the solid lines. The evolutionary tracks in this figure correspond to some of the blue tracks 
in Figure 7, and they all start with a zero-age MS star with mass 1.6 Mo orbiting a point-mass WD with 
mass 0.53 Mo, but each one of them start at a different orbital period. Before the onset of CE evolution, the 
atmospheric Roche lobe overflow model provides accretion rates that can only explain SySts that are close 
filling their Roche lobe, i.e. Ra z 0.90 RRL. 

Moreover, the SySt lifetime predicted by this model is ~ 0.01 — 0.1 Myr. This is orders of magnitude 
shorter than those typically assumed in studies that try to estimate the expected number of SySts in the Milky 
Way [e.g. 333]. In general, regardless of the assumptions, these investigations predict much more Milky Way 
SySts than observed. This problem could be partially/completely overcome with the atmospheric Roche lobe 
overflow model, if the intrinsic population of S-type SySts in the Milky Way host red giants close to Roche 
lobe filling. However, if only a small fraction in the intrinsic population shows evidence for ellipsoidal 
variability, like currently observed [331], then this model most likely does not correspond to the dominant 
mode of mass transfer for these systems. In what follows we discuss the case in which the red giant is 
strongly under-filling its Roche lobe. 


6.3.2 Gravitationally focused wind accretion 


In the case of low-mass red giants, the mass transfer rates provided by Equation | in the atmospheric 
Roche lobe overflow model are only relevant if Ra = 0.85 — 0.9 RRL. If this condition is not satisfied, then 
the accretion rates onto the WD are negligible, and in turn this model cannot explain the observed WD 
luminosities in SySts. This implies that mass transfer must be associated with accretion of part of the stellar 
winds from the red giant, i.e. 


Ma = —T)wind Mwind > (27) 


where wina is the orbit-averaged wind accretion rate efficiency. Observations suggest that the wind mass 
loss rate from a typical red giant donor in SySts is = 1077 Mo yr™! [e.g. 334], implying that mwina Z 0.1 to 
reach the accretion rates of = 1078 — 1077 Mo yr™! needed to explain the hydrogen nuclear burning. 

The most simple case of wind accretion occurs when we assume that stellar winds are supersonic and 
spherically symmetric, which is the so-called Bondi-Hoyle-Lyttleton model [335, 336, 337]. In this case, the 
WD can accrete part of the stellar wind as it orbit through it and the accretion rate efficiency, i.e. the fraction 
of the stellar winds that is effectively accreted by the WD, can be estimated as 


ou. Ý 21-3/2 
Qwind a Vorb 
indBHL = 1 +{—— <l, (28) 
sin 2V1 — e2 (S | (e) | 
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Fig. 9 Predicted accretion rate (M,) onto the WD during pre-SySt and SySt evolution, as a function of the orbital period Por, 
according to the atmospheric Roche lobe overflow model (Equation 1), the Bondi-Hoyle-Lyttleton model (Equations 27 and 28, 
with @wind = 1.5), the wind Roche lobe overflow model (Equations 27, 29 and 30, with Taust = 1500 K and p = 1), and the wind 
compression model (Equation 31, assuming 5 < p°/p%®™ < 15). The grey thick horizontal line marks the minimum accretion 
rate such that they system can be considered a SySt, which is ~ 1078 Mo yr7!. The evolutionary tracks were computed by us 
using the mesa code and correspond to some of those shown in Figure 7. For these selected tracks, we start the simulations 
with a zero-age MS star of mass 1.6 Mo orbiting a point-mass WD with mass 0.53 Mo, which are consistent with the average 
values for the red giant and WD in SySts [21]. We chose nine different initial orbital periods, from left to right: 300, 450, 600, 
800, 1000, 1700, 2400, 3100 and 3800 days, and evolve each one of them until the onset of CE evolution. In all tracks, the 
orbital period increases/decreases as the SySt evolves mainly as a response to red giant donor spin evolution, since it quickly 
becomes synchronized with the orbit while ascending the FGB. For the three binaries with shortest initial orbital period, Roche 
lobe filling occurs when the donor is on the FGB, while for the remaining binaries it takes place when the donor is on the 
AGB. In case the orbital period during the SySt phase is shorter than 10° days, we computed the accretion rates provided by 
the atmospheric Roche lobe overflow model (solid lines) and those obtained with gravitationally focused wind accretion, for 
which we assumed either the Bondi-Hoyle-Lyttleton model, if the donor is on the FGB (dashed blue lines), or the wind Roche 
lobe overflow model, if the donor is on the AGB (dashed red lines), or the wind compression model (green areas). Furthermore, 
we only show the accretion rates provided by the wind Roche lobe overflow model model when the orbital periods are longer 
than 10° days. The average Bondi-Hoyle-Lyttleton and wind Roche lobe overflow accretion rate efficiencies ()windBHL and 
Nwind,wRLOF) are indicated in the figure, while those for the compression wind model are given by (5— 15) X wind BHL- It is also 
indicated the typical SySt lifetime (7), based on the atmospheric Roche lobe overflow and wind compression model, for orbital 
periods shorter than 10° days, and based on the wind Roche lobe overflow model, for orbital periods longer than that. SySts with 
orbital periods shorter than 10° days cannot be easily explained by the Bondi-Hoyle-Lyttleton model, nor the wind Roche lobe 
overflow model, as these mass transfer modes do not provide accretion rates high enough to trigger the symbiotic phenomenon. 
However, the required rates in this period range can be obtained with the wind compression model, if the accretion efficiency 
is 2 5 times higher than that from the Bondi-Hoyle-Lyttleton model. On the other hand, the wind Roche lobe overflow model 
is clearly able to explain SySts with AGB donors at orbital periods longer than 10° days. Finally, the atmospheric Roche lobe 
overflow model always provide sufficiently high accretion rates, when Ra 2 0.85 — 0.9 Ret. 
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where Vwind and vorb are the wind and orbital velocities, respectively, and Qwing is a parameter to be adjusted 
based on observations. 

Depending on the properties of the red giant and the orbital separation, the winds can be gravitationally 
focused towards the orbital plane and more easily accreted by the WD. Even though there are a few models 
available to estimate gravitationally focused wind accretion rates [e.g. 172, 174, 338], we will only consider 
two models to illustrate how this mechanism can lead to accretion efficiencies much higher than those ob- 
tained with the Bondi-Hoyle-Lyttleton model. These two models are (1) the wind Roche lobe overflow model 
[170, 171, 172] and (ii) the wind compression model [195, 174], and they look very promising in the context 
of SySts. 

Let us start with the wind Roche lobe overflow model. During the AGB phase, winds are most likely 
driven by a combination of pulsation-induced shock waves and radiation pressure on dust grains [339]. 
Stellar pulsation and convection induce strong shock waves in the extended atmosphere, pushing the gas 
outwards, which reaches cooler regions. Since most of the gas in a cool AGB star is in molecular form, some 
species may condense into dust grains, giving rise to a dusty shell around the star. These dust grains are then 
accelerated outwards by radiation pressure and the gas is dragged along. In the presence of a companion, the 
gas can remain within the AGB Roche lobe, if the wind acceleration radius is larger than the AGB Roche 
lobe radius [172]. When this happens, instead of having the red giant effectively filling its Roche lobe due 
to nuclear evolutionary processes, the slow winds fill its Roche lobe, enhancing in turn the wind accretion 
onto the WD [170, 171, 173, 175, 338]. This implies that the accretion rates predicted by the wind Roche 
lobe overflow model are expected to be higher than those provided by the Bondi-Hoyle-Lyttleton model, but 
smaller than those obtained with the atmospheric Roche lobe overflow model. 

The accretion rate efficiency in the wind Roche lobe overflow model can be expressed as [172] 


2 
25 R R 
Twind,wRLOF = e (2)|-o284( au +0918( ast) 0.234 <0.5, (29) 


RRL RRL 


where the maximum value 0.5 comes from the results of hydro-dynamical simulations [170, 171], and Rgust 
is the dust condensation radius, corresponding the radius of the wind acceleration zone, given by [340] 


-(4 2 
fa) Taust (4+p)/ 
2 Tet 


Ras = ( (30) 


where Taust is the dust condensation temperature and p is a parameter characterising wavelength dependence 
of the dust opacity. 

Another way to enhance wind accretion is by means of asymmetric radiation-driven winds from rotating 
red giants. Based on the wind compression disk model [195], which explains the equatorial circumstellar 
disks in rapidly rotating Be stars, it has been proposed that gravitational focusing in SySts can be achieved 
if the radiation-driven wind from the red giant donor is not spherically symmetric and is partially confined 
to the equatorial plane around the donor [174]. Compared to Be stars, normal red giants are slow rotators. 
Despite that, red giants in SySts rotate faster than isolated stars with comparable spectral type [341, 342], 
and this mechanism can still work, leading to the formation of a circumstellar gaseous disk-like component 
around the red giant donor. Therefore, part of the wind leaving the red giant can be compressed to its 
equatorial plane, which increases the mass loss around the equatorial plane, proving then accretion rate 
efficiencies that are enhanced in comparison with the spherically symmetric case, i.e. the Bondi-Hoyle- 
Lyttleton model. 
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The required assumption that winds are asymmetric seems to be consistent with observations. For in- 
stance, based on the velocity profile of the wind from the giant near the orbital plane of the eclipsing SySts 
EG And and SY Mus, it was possible to conclude that their winds are not spherically symmetric [343]. In 
addition, the properties of the nebular [O m] 25007 line in EG And provides further support for the wind 
being focusing towards the orbital plane [344]. 

The accretion rate efficiency in the wind compression model can be expressed as [174] 


p(r) ) 


pPr(r) a 


TwindWC = nina ( 


where r is the distance between the center of the red giant and the WD, p‘?"(r) is the density of the spherically 
symmetric wind at a distance r, and p‘(r) is the local density of the compressed wind at a distance r in the 
equatorial plane. The quantity p°(r)/p*?'(r) also depends on the radius and rotational velocity of the red 
giant, on the initial and terminal wind velocities, and on how the wind is accelerated. For typical rotational 
velocities of red giants in SySts (~ 6-10 km s7! [342]) and wind terminal velocities (~ 20 — 50 km s~!), the 
accretion rate efficiencies obtained with the wind compressional model, in comparison with those from the 
Bondi-Hoyle-Lyttleton model, can be enhanced by a factor of ~ 5— 15 [174]. 

We compare in Figure 9 the accretion rates onto the WDs in SySts predicted by the atmospheric Roche 
lobe overflow model, the Bondi-Hoyle-Lyttleton model, the wind Roche lobe overflow model, and the wind 
compression model. It is clear from the figure that the Bondi-Hoyle-Lyttleton model cannot account for 
the accretion rates in SySts. On the other hand, the wind Roche lobe overflow model can reasonably well 
explain the SySts with orbital periods longer than a few thousand days. However, for those systems with 
orbital periods shorter than ~ 10° days, this model is most likely unable to provide the required accretion 
rates, since in this orbital period range it provides accretion rate efficiencies comparable to the Bondi-Hoyle- 
Lyttleton model [172]. The wind compression model is a good candidate to explain SySts in this range, as 
it provides accretion rates = 1078 Mo yr7!. As final comment, the SySt lifetimes predicted by these gravita- 
tionally focused wind accretion models (~ 0.1 — 1 Myr) may also help to bring into agreement predicted and 
observed numbers of Milky Way SySts, although further and more detailed calculations are needed to better 
address this issue. 


6.4 Supergiant and Wolf-Rayet high-mass X-ray binaries 


We have just describe the modes of mass transfer in SySts and SyXBs. Given that sg-HMXBs and WR- 
HMXBs are the massive counterparts of these systems, the same modes could be in principle directly appli- 
cable to sg-HMXBs and WR-HMXBs. 

Let us start with ss-HMXBs hosting NSs in which the supergiants are nearly filling their Roche lobes, 
such as Cen X-3, SMC X-1, LMC X-4 [253], and XMMU J013236.7+303228 [257]. The mass transfer 
in these systems can be described by the atmospheric Roche lobe overflow mode in which the extended 
photosphere of the supergiant is in contact with the Roche radius. In this case, matter flows to the compact 
object through the inner Lagrangian point and an accretion disk around the NS is formed. Given that the 
donor masses in these systems are several times higher than the masses of their NS companions, mass transfer 
will quickly become dynamically unstable mass transfer, triggering CE evolution. The most likely outcome 
of this situation is a complete coalescence of the system, leading to the formation of a Thorne-Zytkow 
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object [345], which would be a very cool red supergiant with a degenerate neutron core, although these 
objects remain thus far hypothetical. 

In other sg-HMXBs and WR-HMXBs, their donors are under-filling their Roche lobes, and the most 
likely mode of mass transfer for them is gravitationally focused wind accretion. An important implication of 
this mass transfer mode is that it can naturally explain wind-fed ultraluminous X-ray sources (ULXs [346]), 
such as the WR-HMXB M101 ULX-1 [202]. In addition, it also offers an explanation for the existence of 
accretion disks in the WR-HMXBs IC 10 X—1 and NGC 300 X-1, since this model is consistent with the 
formation of a wind-captured disk [347]. Moreover, results from recent binary models indicates that the 
contribution of wind-fed ULXs to the overall ULX population is up to ~ 75 — 96% for young (< 100 Myr) 
star-forming environments and up to ~ 49 — 87% for environments like the Milky Way disk [348]. All this 
implies that Roche lobe filling is not the only mode of mass transfer able to provide large accretion rates 
onto the NS/BH in ULXs. 

As a final comment, there is potential important application of gravitationally focused wind accretion in 
the context of HMXBs, which is related to the formation of Be-HMXBs hosting BHs. Having the uncertain- 
ties in the stability criterion for dynamical mass transfer in mind [201, 125], binary models suggest that these 
systems should form through an episode of CE evolution [200]. In case there is no efficient mass and angular 
momentum transfer from the BH progenitor to the B star, the post-CE binary will host a Wolf-Rayet star 
orbiting a slowly rotating B star, which will not exhibit the typical characteristics of Be-HMXBs. However, 
this problem may be overcome in case of efficient wind accretion prior the CE evolution, such that the B 
star would enter this event already as a rapidly rotating star. Further binary models, with consistent spin up 
prescription, are still needed to better evaluate this possibility. 


7 Problems to Overcome and Possible Solutions 


Multi-wavelength observations of several types of binaries as well as large systematic surveys, such as 
the Sloan Digital Sky Survey, have provided new observational constraints on theories of the formation and 
evolution of accreting compact objects. Together with improved numerical modelling, these constraints have 
led to substantial progress in our understanding of the formation of these fascinating objects which have 
been outlined in the previews sections of this review. Despite these improvements, there are still important 
problems to be solved. In what follows we will list the most critical open questions and potential avenues 
towards obtaining answers. 

For the formation of accreting compact objects, mass transfer from a giant star to its companion star is 
fundamental. We clearly need to better understand the conditions that lead to dynamically stable and unsta- 
ble mass transfer. This will only be possible if we find ways to test which assumptions for the mass transfer 
equation, the radius of the donor and the Lagrangian radius provide more realistic results. Comparison of 
detailed model predictions with observed samples of stars that evolve through dynamically stable mass trans- 
fer or common-envelope evolution may perhaps allow us to reach that goal. Post-common-envelope binaries 
consisting of WDs and MS star companions may be numerous enough and easy enough to characterize to 
provide the required observational constraints. The increasing number of merging binary black holes/neutron 
stars detected by the LIGO and Virgo interferometers as well as combined gravitational wave and elec- 
tromagnetic observations of compact binary mergers will also provide useful and strong multi-messenger 
constraints in the high-mass regime. 

Very much related to the point above, we need to better understand common-envelope evolution. If we 
ever want to understand the relative numbers of accreting compact objects and different binaries hosting 
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compact objects, we need to be able to reliably predict the outcome of common-envelope evolution for a 
given system, which apparently is currently not the case. To reach this goal, on one hand we need to aim at 
developing modelling tools that cover the different physical processes and time-scales involved in common 
envelope evolution. On the other hand, observational surveys that provide clean and ideally volume limited 
samples of close binaries can provide key information for theoretically unconstrained parameters. First steps 
in this direction have been performed for white dwarf binaries with low mass star companions, but large, 
complete, and unbiased samples are still not available. 

Dynamically stable non-conservative mass transfer is potentially the dominant (but most likely not the 
only) formation mechanism for systems hosting white dwarfs with orbital periods between ~ 100 and 
~ 1000 days. These systems are often found in eccentric binaries, which should not be the case if our under- 
standing of tidal evolution is correct. Despite the several proposed mechanisms to enhance the eccentricity 
during dynamically stable non-conservative mass transfer, they all fail to some extent. Either the theory of 
tidal interaction and/or the eccentricity pumping mechanisms need serious revision, or another more promis- 
ing, still unknown, mechanism is needed. On the top of that, we have virtually no idea how much mass and 
angular momentum is lost during dynamically stable non-conservative mass transfer, which makes predic- 
tions rather uncertain. Given that mass loss is a huge problem already in single star evolution, it is unlikely 
that we will soon find the answer by performing theoretical simulations. Decisive progress therefore again 
depends to a significant degree on observational constraints from samples of post-mass-transfer systems. 

A problem equally serious as those listed above is that we do not properly understand orbital angular 
momentum loss in close binaries with both accreting and non-accreting compact objects. Problems related 
very likely to our ignorance of magnetic braking and consequential angular momentum loss are numerous 
when we confront standard prescriptions with observations, e.g. the average white dwarf mass in cataclysmic 
variables is much larger than expected, mass transfer rates in cataclysmic variables and low-mass X-ray 
binaries inferred from observations are very different from those predicted by our models, the predicted 
donor masses in black hole low-mass X-ray binaries are significantly higher than observed, the observed 
orbital period minimum in cataclysmic variables is significantly longer than predicted, the predicted close 
millisecond pulsar binaries come from a range of the parameter space so narrow that their formation becomes 
compromised. 

Promising revisions in our understanding of the orbital angular momentum loss due to magnetic brak- 
ing have been proposed in the past years, such as the convection and rotation boosted prescription, which 
provides stronger orbital angular momentum loss rates and can potentially solve the major problems related 
to the mass transfer rates in low-mass X-ray binaries and consequently the formation of their millisecond 
pulsar binary descendants. In addition, given that a dynamical origin of black hole low-mass X-ray binaries 
appears unlikely, this new prescription perhaps in combination with consequential consequential angular 
momentum loss, e.g. due to circumbinary disks, seems a fruitful path to explore in binary models aiming at 
reproducing their characteristics. However, this prescription cannot be the solution for cataclysmic variables, 
since observations suggest lower orbital angular momentum loss rates than those predicted by the standard 
prescription, which implies a weaker magnetic braking. 

Any important question arising when all these systems are taken into account is whether there is a uni- 
versal magnetic braking prescription able to simultaneously explain all of them. Given the intrinsic different 
donor types in cataclysmic variables and low-mass X-ray binaries, a possible universal prescription could be 
achieved if the donor structure is also taken into account in the formulation of magnetic braking. The over- 
whelming majority of donors in cataclysmic variables are M-type MS stars, while most donors in low-mass 
X-ray binaries seem to be MS stars of earlier spectral types or subgiants. If a magnetic braking prescription 
derived from first principles could provide weaker orbital angular momentum loss rates for MS stars with 
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initial masses smaller than ~ 0.8 Mo and stronger for MS stars initially more massive than that, then such a 
recipe would be a reasonable candidate for a universal magnetic braking prescription. 

All that said, despite the significant progress achieved during decades of investigation, formation and 
evolution of accreting compact objects is a remarkably interesting wide open field. Future observational and 
theoretical efforts will hopefully improve our physical understanding of these systems. 
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